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Abstract 



The baryon asymmetry of the universe can be explained by the 
out-of-equilibrium decays of heavy right-handed neutrinos. We analyse 
this mechanism in the framework of a supersymmetric extension of the 
Standard Model and show that lepton number violating scatterings are 
indispensable for baryogenesis, even though they may wash-out a gener- 
ated asymmetry. By assuming a similar pattern of mixings and masses 
for neutrinos and up- type quarks, as suggested by SO (10) unification, 
we can generate the observed baryon asymmetry without any fine 
tuning, if (B — L) is broken at the unification scale A GUT ~ 10 16 GeV 
and, if m Vfl ~ 3 • 1CT 3 eV as preferred by the MSW solution to the solar 
neutrino deficit. 



Zusammenfassung 

Die Baryonasymmetrie des Universums kann durch den Zerfall schw- 
erer rechtshandiger Neutrinos aufierhalb des thermischen Gleichgewichts 
erklart werden. Wir untersuchen dies im Rahmen einer supersym- 
metrischen Erweiterung des Standard-Modells und zeigen, daB lepton- 
zahlverletzende Streuprozesse, die eine erzeugte Asymmetrie wieder ver- 
nichten konnen, fur die Baryogenese unverzichtbar sind. Nimmt man 
fur Quarks und Leptonen ahnliche Massen und Mischungswinkel an - 
wie von SO(10)-vereinheitlichten Modellen nahegelegt — so kann man 
die beobachtete Baryonasymmetrie erzeugen. Dazu wird (B — L) an 
der Vereinheitlichungsskala A GUT ~ 10 16 GeV gebrochen, und m„ ~ 
3 • 10~ 3 eV angenommen, wie es die MSW-L6sung des solaren Neutrino- 
Problems nahelegt. 
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Introduction 



The observed baryon asymmetry of the universe is one of the most intriguing problems of particle 
physics and cosmology. This asymmetry, which is usually expressed as ratio of the baryon density 
tib to the entropy density s of the universe, 

y B = M = (0.6 - 1) • 1(T 10 , 

s 

could in principle be an initial condition of the cosmological evolution. However, this is not com- 
patible with an inflationary phase which seems to be required in a consistent cosmological model 
Hence, the baryon asymmetry has to be generated dynamically during the evolution of the 
universe. This is possible if baryon number is not conserved, if C and CP are violated, and if the 
universe is not in thermal equilibrium [0]. 

Although the Standard Model (SM) contains all the necessary ingredients, it is not possible to 
explain the baryon asymmetry within the SM, i.e. one has to envisage extended theories. Grand 
unified theories (GUTs) are attractive for various reasons and there have been many attempts 
to generate Yb at the GUT scale 0. However, these mechanisms are difficult to reconcile with 
inflationary scenarios which require reheating temperatures well below the GUT scale. 

Preheating, i.e. the non-thermal decay of the oscillating inflaton at the end of inflation via para- 
metric resonance @, may re-open the window for GUT baryogenesis, since it enables the coherent 
decay of the inflaton condensate into particles that are more massive than the inflaton itself. How- 
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ever, recent calculations indicate that parametric resonance may be ineffective in most inflationary 
models, if the back reaction of the produced particles onto the condensate, the rescattering of the 
decay products, and the expansion of the universe are taken into account ||. 

In supersymmetric theories, the influence of baryon number carrying scalar condensates along 
flat directions of the scalar potential, i.e. Affleck-Dine baryogenesis ||, requires further studies, 
since it is not clear under which conditions this mechanism can generate a baryon asymmetry of 
the requested magnitude @. 

During the evolution of the early universe, the electroweak phase transition is the last opportu- 
nity to generate a baryon asymmetry without being in conflict with the strong experimental bounds 
on baryon number violation at low energies ||. However, the thermodynamics of this transition 
indicates that such scenarios are rather unlikely ||. 

Therefore, the baryon asymmetry has to be generated between the reheating scale and the 
electroweak scale, where baryon plus lepton number (B + L) violating anomalous processes are in 
thermal equilibrium thereby making a (B — L) violation necessary for baryogenesis. Hence, 
no asymmetry can be generated within GUT scenarios based on the gauge group SU(5), where 
{B — L) is a conserved quantity. 

Gauge groups containing SO(10) predict the existence of right-handed neutrinos. In such the- 
ories (B — L) is spontaneously broken, one consequence being that the right-handed neutrinos can 
acquire large Majorana masses, thereby explaining the smallness of the light neutrino masses via 



the see-saw mechanism [11]. Heavy right-handed Majorana neutrinos violate lepton number in 
their decays, thus implementing the required {B — L) breaking as lepton number violation. This 
leptogenesis mechanism was first suggested by Fukugita and Yanagida and has subsequently 



been studied by several authors (see, e.g., refs. [13-21]]). As detailed studies have shown, the ob- 
served baryon asymmetry can be generated in non-supersymmetric [0, 14,ll| and supersymmetric 
theories [lq.17]. 

If one assumes a similar pattern of mass ratios and mixings for leptons and quarks and, if 
m u ~ 3 • 10 -3 eV as preferred by the MSW solution to the solar neutrino problem, leptogenesis 



implies that (B — L) is broken at the unification scale [15]. This suggests a grand unified theory 
based on the group SO(10), or one of its extensions, which is directly broken into the standard 
model gauge group at the unification scale ~ 10 16 GeV. However, for a successful gauge coupling 
unification, such a GUT scenario requires low-energy supersymmetry. 



Supersymmetric leptogenesis has already been considered in refs. [16,18] in the approximation 
that there are no lepton number violating scatterings which can inhibit the generation of a lepton 
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number. Another usually neglected problem of leptogenesis scenarios is the necessary production 
of the right-handed neutrinos after reheating. In non-supersymmetric scenarios one has to assume 
additional interactions of the right-handed neutrinos for successful leptogenesis jjlJJ] . 

In this thesis, we investigate super symmetric leptogenesis within the framework of the mini- 
mal supersymmetric standard model (MSSM), to which we add right-handed Majorana neutrinos, 
as suggested by SO(10) unification fi~7|| . Since CP asymmetries in the decays of these neutrinos 
are one of the principal ingredients of this model, we start by considering possible sources of CP 
violation in decays of Majorana neutrinos in the next chapter. In addition to the usually considered 
one-loop vertex corrections we show how self-energy contributions to the CP asymmetry, which 
have previously been considered in refs. [18-^1]], can be consistently taken into account [22]. For 
simplicity we only consider the non-supersymmetric leptogenesis scenario. However, our results are 
easily generalized to the supersymmetric case. 

In chapter || we present superfield techniques, which simplify calculations in theories with ex- 
act supersymmetry These techniques are used in chapter ^ where we introduce supersymmetric 
leptogenesis. In particular, we discuss the neutrino decays and scattering processes that one has 



to take into account to be consistent [17]. In chapter [| we develop the full network of Boltzmann 
equations necessary to get a reliable relation between the input parameters and the final baryon 
asymmetry. We work out the parameter dependence of the generated baryon asymmetry, and show 
that by neglecting the lepton number violating scatterings one largely overestimates the generated 
asymmetry, and that in our scenario the Yukawa interactions are strong enough to produce a ther- 
mal population of right-handed neutrinos at high temperatures. Finally, we show in chapter || that 
by assuming a similar pattern of masses and mixings for leptons and quarks one gets the required 
value for the baryon asymmetry without any fine tuning, provided (B — L) is broken at the GUT 
scale, and the Dirac mass scale for the neutrinos is of order of the top quark mass, as suggested by 
SO (10) unification 

In appendix [A| we summarize some standard formulae for one-loop integrals. In appendix |B] we 
introduce our spinor notation and compile formulae which are needed for the superfield calculations 
of chapters Q and |3|, while the Feynman rules for component field calculations are presented in 
appendix |C|. After a brief review of thermodynamics in an expanding universe in appendix [D|, we 
present the cross sections for the scattering processes discussed in chapter |3] in appendix ^. Finally, 
in appendix |F] we discuss some limiting cases in which the corresponding reaction densities can be 
calculated analytically. 



Chapter 1 

CP Asymmetry in Majorana 
Neutrino Decays 



In this chapter we study how self-energy diagrams can be consistently taken into account when 
computing CP asymmetries in heavy particle decays. This is not obvious, since the naive pre- 
scription leads to a well-defined result for the CP asymmetry, whereas the individual partial decay 
widths are infinite. 

We investigate this problem in the case of heavy Majorana neutrinos, which are obtained as mass 
eigenstates if right-handed neutrinos are added to the standard model. Since they are unstable, they 
cannot appear as in- or out-states of S-matrix elements. Rather, their properties are defined by S- 
matrix elements for scatterings of stable particles mediated by the unstable neutrino |23|] . By using a 
resummed propagator for the intermediate neutrino, we can separate two-body scattering processes 
in resonance contributions and remainder. While the CP asymmetries of two-body processes 



vanish [24], the resonance contributions yield a finite CP asymmetry which can be assigned to the 



intermediate neutrino 

1.1 Self-energy and vertex corrections 

We consider the standard model with three additional right-handed neutrinos. The corresponding 
Lagrangian for Yukawa couplings and masses of charged leptons and neutrinos reads 

C Y = kH\te R + keH j \tv R -^Mv R + h.c. , (1.1) 

where Z L = (^ L ,e L ) is the left-handed lepton doublet and H = (H + ,H°) is the standard model 
Higgs doublet. A/, X u and M are 3x3 complex matrices in the case of three generations. One 
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Figure 1.1: Leading order contributions to the self-energy of the heavy Majorana neutrinos. The 
fermion flow has been chosen parallel to the external momentum q. 

can always choose a basis for the fields f R such that the mass matrix M is diagonal and real with 
eigenvalues Mj. The corresponding physical mass eigenstates are then the three Majorana neutrinos 
Ni = + v^L. At tree level the propagator matrix of these Majorana neutrinos reads 

'*<«> = J^wm- (L2) 

This propagator has poles at q 2 = M 2 corresponding to stable particles, whereas the physical 
Majorana neutrinos are unstable. This is taken into account by summing self-energy diagrams in 
the usual way, which leads to the resummed propagator 

> s ^ = 4-M-m - (L3) 



At one- loop level the two diagrams in fig. 1.1 yield the self energy^] 



£l>) = {4P*U S «V) + W^W *W) , (1.4) 

where P R>L = | (1^75) are the projectors on right- and left-handed chiral states. S R and S L are the 
contributions of the diagrams figs. (|l.l| a) and (|l.l| b), respectively. They can be written as products 
of a complex function a(q 2 ) and a hermitian matrix K, 

^{q 2 ) = ^ n {q 2 )) T = a{a 2 )K , K = A£ X u . (1.5) 



a 



(q 2 ) is given by the usual form factor Bo(q 2 , 0, 0) defined in eq. (A. 5), whose finite part reads in 



the MS-scheme, 



a(q 2 ) = T ^(lJ-p-2- l ne(q 2 )) . (1.6) 



1 For the calculations we use the non-supersymmetric subset of the Feynman rules in App. by identifying the 
SM Higgs doublet H with the supersymmetric scalar Higgs doublet H2- 
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Figure 1.2: One-loop corrections to the couplings of heavy Majorana neutrinos Nj to anti-lepton 
Higgs states (a) and lepton Higgs states (b). The fermion flow has been chosen according to the 
external lepton lines. 

For simplicity we will often omit the argument of a in the following, however one should keep in 
mind that a depends on q 2 . 

According to eqs. ( |Q| ) and ( |1.4[ ) the resummed propagator S(q) satisfies 

^f(l-S R (g 2 ))P R + (l-E L (g 2 ))P L 



M 



S(q) 



1 . 



The fermion propagator S(q) consists of four chiral parts 

S(q) = P R S™(q 2 ) + P L S L V) + P^S LR (q 2 ) + P^S^tf) . 



[1.7) 



(1.8) 



Inserting this decomposition into eq. (1.7), and multiplying the resulting equation from the left and 
the right with chiral projectors Pr,l, yields a system of four coupled linear equations for the four 
parts of the propagator. The solution reads 

-,2 i-l 



S RR (q 2 ) 
S LR {q 2 ) 
S LL (q 2 ) 
S RL (q 2 ) 



(l-X h (q 2 ))^(l-X R (q 2 ))-M 



1 

M 



(1 - £ R (^))S R V) , 



(l-S R (g 2 ))|:(l-S L (g 2 ))-M 



(1 - £ L (^))S L V) • 



M 



(1.9) 
(1.10) 
(1.11) 
(1.12) 



As discussed below, the diagonal elements of S(q) have approximately the usual Breit-Wigner form. 



In addition to the self-energy we need the one-loop vertex function. The two expressions for 
the coupling of N to I, (fig. |L2|a) and N to I, H (fig. |0|p) can be written as 



f3a,ab 



(q,p) 



{KMb(q,p)Xl) jiq ^ + {KMc(q,p)\ T v ) jiPtl (P^)^ , (1.13) 



,i 1 
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{XlMb{ q ,p)K) iiq ^ + {XlMc{q,p)K) ijP ^ (P K ^)a/3 ■ (1.14) 



Here b(q,p) and c(q,p) are diagonal matrices whose elements are given by the standard three point 
form factors Co and C\ 2 defined in appendix |A.3 , 

1 



h(q,p) 

Ck(q,p) 



16vr 2 
1 

16vr 2 



[C (-p - q, q, M k , 0, 0) + C 12 (-p - q,q,M k ,0,0)] , 



[C (-p - q, q, M k , 0, 0) + 2C 12 (-p - q, q, M k ,Q, 0)] 



(1.15) 
(1.16) 



Since we shall only consider amplitudes with massless on-shell leptons, the terms proportional to 
c k will not contribute. We shall only need the imaginary part of b k which is given by 



1 / M 2 \ 

where the function / is defined as 

/(x) = V^l-(l + z)m(^^) . 

1.2 Transition matrix elements 



(1.17) 



(1.18) 



The two lepton-number violating and the two lepton-number conserving processes are shown in 
figs. |1.3|a- |1.3| d. Consider first the contributions of the full propagator, where the full vertices are 
replaced by tree couplings. The four scattering amplitudes read 



(W)H e ^q-p')\l?(p)H b (q-p)) 



^(p')H e (q-p')\lt(p)H b ^q-p)) 



(l'*(p')H e (q-p'M(j ) )H b (q-p)) 



(lUp')H e \q-p')\lt(p)H b \q-p)) 



hi 

(CP L v(p')) T St k L (q)(PMp)), (1-19) 

+ie a bede(Xl)ij(Xi)ki 

(u(p')P R )S^(q)(v(p)P R Cf , (1.20) 

(u(p')P R )S%(q)q'(PMp)), (1-21) 

hi 

(CP l v(p')) T SSf (?) 4 (Hp)P«C) t . (1.22) 



Here a, b, d, e denote the SU(2) indices of lepton and Higgs fields and k, I are generation indices. 
The relative signs follow from Fermi statistics. 
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We are particularly interested in the contributions of a single heavy neutrino to the scattering 
amplitudes. In order to determine these contributions we have to find the poles and the residues of 
the propagator matrix. Here an unfamiliar complication arises due to the fact that the self-energy 
matrix is different for left- and right-handed states. Hence, the different chiral projections of the 
propagator matrix are diagonalized by different matrices. 

5 LL and S Rn are symmetric complex matrices, since £ L (<7 2 ) = (S R (g 2 )) T . Hence, S LL and S RR 
can be diagonalized by complex orthogonal matrices V and U, respectively, 

S LL (q 2 ) = V T {q 2 )MD{q 2 )V{q 2 ) , S m {q 2 ) = U T (q 2 )MD(q 2 )U(q 2 ) . (1.23) 

Splitting the self-energy into a diagonal and an off-diagonal part, 

S L (g 2 ) = S D (g 2 ) + S N ( (? 2 ), (1.24) 

one finds 

D~\q 2 ) = 9 2 U - £ D (9 2 )) 2 ~M 2 + 0(S 2 ) . (1.25) 

One can easily identify real and imaginary parts of the propagator poles. The pole masses are given 
by 

M- = Z M M 2 , Z M . = 1 + (in ^- - 2^ , (1.26) 

and the widths are Tj = KaMi/(8ir). In the vicinity of the poles the propagator has the familiar 
Breit-Wigner form 



Zi _ . Ku f. Mf 



q 2 - M\ + iMiTi 8vr 2 V V 2 J 

We can now easily write down the contribution of a single resonance Ni with spin s to the 
lepton-Higgs scattering amplitudes. Suppressing spin indices for massless fermions, one has 



(J^p'^iq-p'Mip^iq-p^ = (J^p'^iq-p'm^s)) 

zD^q 2 ) (N l (q, S )\lt(p)H\q-p)) , (1.28) 

(l$(j/)H*(q-i/)\lf(p)H>*(q-p)) l = {l$(j/)H°(q - j/)|JV,(?, *)> 

iDtiq 2 ) (N^sMip^iq-p)) , (1.29) 

{l d ] {p l )H e (q-p'M{p)H\q-p)) l = (lf( P ')H e (q - j/)|JVj( g , s)) LC 

iDtiq 2 ) (NfasMWH^q-p)), (1.30) 

{lfa)H*(q-i/)\lf(p)H*(q-p)) l = (Jfa'^iq-p'm^s))^ 

iDtiq 2 ) (N l (q, S )\lUp)H b \q-p)) . (1.31) 
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It v ^ 



(a) 



H 



b/ 1 







Figure 1.3: s-channel contributions to lepton-Higgs scattering, including full propagators and ver- 
tices. 

Here the subscript LC distinguishes an amplitude defined by a lepton-number conserving process 
from the same amplitude defined by a lepton-number violating process. From eqs. (tL9l)-([Tl2D and 
(|L2|) one finds 



^{p'W^q-p'm^s)) 
(N^sMip^iq-p)) 



(l^H^q-p'm^s)) 



{l^H^q-p'm^s^Lc 



+ie a b{V{q 2 )\ T u ) H u.(qf,M,)P L u(p), 
-ie de (A„y% 2 ))^ s (g,M z )P L ?;(p') 
-itab (u(q 2 )Xl) . v(p)P K v a {q,Mi) , 



+ie de [KU T {q 2 )) , l u{p')P J!L u a { k q,Mi) 



(W^iq-p'm^s)) 



LC 



+ie de [K-{1 - Z h (q 2 ))V T (q 2 )M 

-ie de - ^{q 2 ))U T {q 2 )M 

Vsiq^M^PMp') , 



jl 



jl 



(1.32) 
(1.33) 
(1.34) 
(1.35) 

(1.36) 
(1.37) 



where we have used the identity CvJ(p) = u s (p). 

Eqs. ( |l.32| ) and ( 1.33 ) describe the coupling of the Majorana field N to the lepton fields Zj and 
the Higgs field H, and eqs. ( 1.34 ) and ( 1.35|) give the couplings of N to the charge conjugated 



fields U and W. In the case of CP conservation, one has \ Vi - = A^., which implies K = K T and 
therefore 



S L (g 2 ) = E R (q 2 ) , V(q 2 ) = U(q 2 ) . 



(1.38) 
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This yields 



(Ni(q,s)\lt(p)H b (q-p)) = (N t (q, s)\l?(p)H b \q - p)) 



(1.39) 



with q = (qo, —q), p = (po, —p), as required by CP invariance. 

The amplitudes given in eqs. (|1 .32 ) - ( 1.35| ) have been obtained from the lepton-number violating 
processes figs. L3a and |1.3| b. The lepton-number conserving processes figs. |1.3| c and |1.3| d yield the 
amplitudes given in eqs. ( 1.36| ) and ( 1.37 ). The consistent definition of an on-shell contribution of a 
single heavy Majorana neutrino to the two-body scattering amplitudes requires that the transition 
amplitudes extracted from lepton-number conserving and lepton-number violating processes are 
consistent. This implies 



(J^p'^iq-p'm^s)) = (J^p'^iq-p'm^sVLc. 



(1.40) 
(1.41) 



From eqs. Q1.33D , ( |1.35| ), ( |1.36| ) and (1.37) it is clear that these relations are fulfilled if the mixing 
matrices V(q 2 ) and U(q 2 ) satisfy certain consistency relations. Assuming that the matrix X u has 
an inverse, one reads off 



Uij(M?) = (MV(Mi)(l-Z n (Mf))± 



Vij(M?) = (mXJ{M?)[1-^{M1))^ 



(1.42) 



(1.43) 



The matrices V and U are determined by the requirement that the expressions (cf. eqs. (|L23|) ) 



V(q 2 ) (^V))"* V T (q 2 ) = V(q 2 ) ((l - S^ 2 )) ^ (l - £ L (g 2 )) - m\ V T (q' 

Utf) (^V))" 1 U T (q 2 ) = U(q 2 ) ((1 - S L (, 2 )) g (1 - £ R (g 2 )) - M) U T {q 
are diagonal on-shell, i.e., at q 2 = M 2 . Using S L = S D + S N , and writing 

V(q 2 ) = l + v(q 2 ), v(q 2 ) = -v T (q 2 ), 
U(q 2 ) = l + u(q 2 ), u{q 2 ) = -u T {q 2 ), 
a straightforward calculation yields 

v i3 (q 2 ) = ^(q 2 ) (M i S N ..(( Z 2 )+M j S N . j (g 2 )) , 
Uijiq 2 ) = Wlj (q 2 ) (M^iq 2 ) + M^ Nji {q 2 



(1.44) 
(1.45) 

(1.46) 
(1.47) 

(1.48) 
(1.49) 
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where 



Wij(q 2 



(Mi - Mj 



MiMi 
1 + -Ar 1 - 



2a(q 2 ) [MiK 



j j 



MjK u ) 



(1.50) 



These equations give the matrices V and U to leading order in E N . They are meaningful as long 
as the matrix elements of E N are small compared to those of to . 

Inserting eqs. (|L48p and (|L"49|) m eqs. flOp and ( [Ug ), one finds that the consistency conditions 
for the mixing matrices V and U are fulfilled to leading order in E N . We conclude that the 
contribution of a single heavy neutrino to two-body scattering processes can indeed be consistently 



defined. The pole masses are given by eq. (1.26) and the couplings to lepton-Higgs initial and final 
states are given by eqs. (1.32D-(fL~35). 



1.3 CP asymmetry in heavy neutrino decays 



It is now straightforward to evaluate the CP asymmetry in the decay of a heavy Majorana neutrino, 

T(Ni -> IH) - T(Ni -» Lfft) 



1 F(Ni IH) + T(Ni -► [Eft) • 
From eqs. (|L3|) and (|L3|) one obtains for the partial decay widths, including mixing effects 



(1.51) 



2 



T M (Ni^W) oc ^|(A^ T (Mf)) 

3 

T^Ni^lH) cc £|(A*17 T (M?)) 



To leading order in X 2 this yields the asymmetry (cf. eqs. (|1.46|) , (|1.47| )), 

ef = l-ReiiuiMfWu - (v(Mf)K T ) u } . 



(1.52) 
(1.53) 



Using eqs. flPsj ) - ( |L50| ) and (|E|), one finally obtains 



(1.54) 



(1.55) 



Consider first the case where differences between heavy neutrino masses are large, i.e., |Mj 



Mj\ 3> |r, — Tj\. Eq. ( |1.55[ ) then simplifies to 



-V 

»7T ^ 



MiMj M^n -} 



(1.56) 



This is the familiar CP asymmetry due to flavour mixing fig] . It has previously been obtained 
by considering directly the self-energy correction to the Majorana neutrino decay, without any 
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Figure 1.4: u-channel contributions to lepton-Higgs scattering. 

resummation. The CP asymmetry £j reaches its maximum for |Mj — MA ~ |3?^ — Tj\, where the 
perturbative expansion breaks down. 

Interesting is also the limiting case where the heavy neutrinos become mass degenerate. From 
eq. ( |1.55| ) it is obvious that the CP asymmetry vanishes in this limit. The vanishing of the CP 
asymmetry for mass degenerate heavy neutrinos is expected on general grounds, since in this case 
the CP violating phases of the matrix K can be eliminated by a change of basis. 

The CP asymmetry due to the vertex corrections is easily obtained using eqs. ( 1.13|) , ( |l,14j ), 
(|1.17|) and ( |1.18 ). The partial decay widths corresponding to the full vertex read 



T v (Ni^lH^) oc J2\(X U (1- MbK T M)) 

3 

T v {Ni^lH) oc ^2\{X* U (1- MbKM)) 



2 

3%\ i 



For the corresponding CP asymmetry ( |1.51 ) one obtains the familiar result 

'2 1 /* f 2" 



I _Im{^.J 
8^r 



Stt K :: J \ M 2 



where the function f(x) has been defined in eq. ( |1.18 ) 



(1.57) 
(1.58) 



(1.59) 



1.4 CP asymmetries in two-body processes 

Let us now consider the CP asymmetries in two-body processes. Here we have to take into account 
the s-channel amplitudes shown in figs. L3a and |l.3| b, with vertex functions up to one- loop, and 



the two u-channel amplitudes depicted in figs. |1.4j and L4 b. For the u-channel amplitudes vertex 
and self-energy corrections can be omitted to leading order since the absorptive parts vanish. 
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In the following we shall evaluate various contributions to the CP asymmetry 

= A|.M| 2 = \M(IH^IH)\ 2 -\M(IH^IH^)\ 2 
£ 2\M\ 2 \M(lHi->lH)\2 + \M(lH -»rfft)|a ' l ' ' 

where we always sum over generations in initial and final states. There are contributions from 
the full s-channel propagator, A|A^|g, from the interference between s-channel amplitudes at tree- 
level and with one-loop vertex corrections, A|.M| 2 r , the interference between tree- level s-channel 
and u-channel amplitudes, A|.M| 2 n , and the interference between s-channel with one- loop vertex 
corrections and u-channel amplitudes, A|.M| 2 r . 

Consider first the CP asymmetry e s due to the full propagator. The contribution of a single 
intermediate neutrino JVi is (cf. Q, flD% ), ( |L33| )) 



\M t (lH^lH^\^\D t ( q 2 )\ 2 Y,\(y(l 2 )^h\ 2 El( A - yT ^ 2 ))^| 2 - (1-61) 

j k 

Comparison with eq. ( |1.52| ) yields immediately 

\Mi(lH -> Lfft)|2 K | jD .( 9 2)|2 rM ( iV . ^ ^t)2 _ (L62) 

Similarly, one has for the charge conjugated process 

\Mi(JHi - IH)\ 2 S oc |A(9 2 )| 2 r M (iVi - ZF) 2 . (1.63) 

The corresponding CP asymmetry is, as expected, twice the asymmetry in the decay due to mixing, 

(i) _ k\Mj\ 2 s M 
£s ~ 2\M t \ 2 ~ 1 ■ 1 j 

It is very instructive to compare the contribution of a single resonance with the CP asymmetry 
e s for the full propagator. Due to the structure of the propagators S LL and S RR it is difficult to 
evaluate e s exactly. However, one may easily calculate e s perturbatively in powers of S N , like the 
mixing matrices V(q 2 ) and U{q 2 ) in the previous section. 

The full propagator (cf. ( |1.7|) ) reads to first order in S N , 



S(q) = S D (q) + S D (q)<i [S^(g 2 )P R + S N (g 2 )P L ] S D (q) + ■■■ , (1-65) 

where (cf. JL25D ) 

S D (q) = W ~ ^(q 2 )) + Ml D(q 2 ) . (1.66) 
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It is now straightforward to calculate the matrix elements of the two-body processes, summed 
over generations in initial and final states, 

\M{IH^W)\ 2 S = Wp-p'q 2 ! -^Tr [KM D(q 2 )K T MD* {q 2 )] + 

\2q 

Re{Tr[KMD(q 2 )Zl(q 2 )(l - Z D (q 2 ))D(q 2 )K T MD*(q 2 ) + 

KD(q 2 )(l - Z D (q 2 ))^(q 2 )MD(q 2 )K T MD*(q 2 )\ } + . . . J . (1.67) 

This yields for the sum and the difference of the CP conjugated matrix elements, 

2\M\ 2 = 16p-p'Y,A ij + ... , (1.68) 

y 

A\M\ 2 = -16p-p l Y / (B ij + C ij ) + ... , (1.69) 



where 



A i3 = ReiKfjMiMjDjiq^D^q 2 )} , (1.70) 
By = ilmlK^MiMjDjiq^DUq 2 ) , (1.71) 
C tj = 4q 2 Re{ia(q 2 )Im{K 2 ..}M i M j (l - S^ 2 ),) K^(g 2 )|A(g 2 )| 2 } . (1.72) 
For q 2 ~ Mf the expressions Au and Cy are dominated by the contribution of a single resonance 



A ti ~ K&MftDitf)?, (1.73) 
- i l ^ K ^^M-2Ku\DM 2 )? ■ (1-74) 

From eqs. (1.56), ( 1.73j ) and ( |1.74| ) one reads off that the sum over Cjj yields precisely the contri- 
bution of the resonance iV, to the CP asymmetry, 

_ £i ^ _ _ J_ v ^ Im{ ^ } - 2,m n 7 ^ 

The second contribution to the CP asymmetry e s is due to the sum over Bij (cf. eq. ( 1.71D ). 
By involves two different propagators (i ^ j) and corresponds to an interference term. Using 
D*- 1 {q 2 ) = q 2 - M 2 - 2a*(q 2 )q 2 K jj and 2q 2 lm{a(q 2 )}K ii = -Jm{DT x {<?)} , one can rewrite dj 
as follows, 

Ca = -i1m{K* JMiMjD^q^Dr^q^Diiq^Djiq 2 )] 2 . (1.76) 
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Comparing eqs. ( 1.71] ) and ( 1.76| ) it is obvious that the sum of both terms, i.e., the CP asymmetry 



e s corresponding to the full propagator, is identically zero! The pole contribution is cancelled by 
the interference of the pole term with an off-shell propagator. 

The contribution to the CP asymmetry A|A4|gp can be computed in a similar manner. The 
diagrams fig. |1.3|a and |l,3|b yield two contributions for the two vertices. After some algebra one 



obtains the result (cf. ( 1.15 )) 



A|M|* r = -Mp-p 1 q 2 Y,D i]k + ... , (1.77) 
D ijk = ImlK^KjkK^Imihiq^MkMjD^D^q 2 ) . (1.78) 

For q 2 ~ M 2 , one reads off that the sum over Dy k yields, as expected, twice the vertex CP 
asymmetry, 

£sHM * ] * —K— f \Ml) ~ ^ ■ (L79) 



A result very similar to eqs. ( 1.77D , ( P--78 ) is obtained for the asymmetry A|A^|^ U , the interference 
between tree-level s-channel and u-channel amplitudes. One finds (u = (q — p — p') 2 ), 

&\M\ 2 S>U = -32 p ■ p' q 2 Y J E ijk + ... , (1.80) 

i,j,k 

E ljk = ImlKikK^K^Imlaiq^MkMjDiiq^D^q^DUu) . (1.81) 
Integrating the expressions over phase space and using 

r^w^^M), (1 . 82) 

J_ q 2 u-Ml M k V 9 2 J 



one finds the cancellation 

f0 



/ du(A\M\l r + A\M\l u ) = 0. (1.83) 

J-q 2 

Finally, we have to consider the CP asymmetry A|A^| 2 r . A straightforward calculation yields 

A|AC,r = -Z2p-p' q 2 Y j F ijk + ... , (1.84) 
F ijk = ImlK.kKjkKjijlmlbkiq^MkMjDiiq^D^u) . (1.85) 

After integration over u the resulting matrix Fij k is antisymmetric in the indices j and k. As a 
consequence, the asymmetry Z\|^W |^ r is identically zero. 
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As we have seen, the total CP asymmetry vanishes to leading order in X 2 . This result has 



previously been obtained in [24|. It follows from unitarity and CPT invariance. The considered 



T-matrix elements satisfy the unitarity relation 

2lm{lH\T\lH) = (IH\T ] T\IH) . (1.86) 

If, in perturbation theory, the leading contribution to the right-hand side is given by two-particle 
intermediate states, one has 



^2(IH\T*T\IH) = {WH\T\IH)\ 2 + \{V H^\T\IH)\ 2 \ + ... . (1.87) 



l l,v 
CPT invariance implies 



(l'H\T\lH) = (lrf\T\l>rf) . (1.88) 
From eqs. ( 1.86j ) - ( |l.88| ) one then obtains 

^ (\{Vrf\T\lH)\ 2 - \{l' H\T\IH^)\ 2 \ + . . . = . (1.89) 



In [24 1 it was concluded that away from resonance poles, where ordinary perturbation theory holds, 



the CP asymmetry ( 1.89 ) vanishes to order A^. Corrections due to four-particle intermediate 
states are O(A^). In this paper we have developed a resummed perturbative expansion in powers 
of S N which is also valid for s ~ M 2 . The same argument then implies that in this case the CP 



asymmetry ( 1.89 ) vanishes to order X 2 whith corrections O(A^). 

The nature of the cancelation is different for different subprocesses. For the full propagator, the 
CP asymmetry vanishes identically for fixed external momenta. Interference contributions between 
various s-channel and u-channel amplitudes cancel after phase space integration. In applications 
at finite temperature the standard practice Q is to treat in the Boltzmann equations resonance 
contributions and the remaining two-body cross sections differently. This procedure yields for the 
CP asymmetry of the decaying heavy neutrino iVj the sum of mixing and vertex contribution, 



Chapter 2 



Perturbation Theory in Superspace 



In this chapter we give a short review of the superspace formulation of theories with global su- 
per symmetry. Since in the following chapters we will only encounter chiral superfields, we restrict 
ourselves to chiral superfields and leave aside the superspace formulation of gauge theories, which 
can be found e.g. in ref. |25|]. In particular, we will present the powerful calculational tool of per- 
turbation theory in superspace, which drastically simplifies calculations of 5-matrix elements in 
theories with exact super symmetry. 



2.1 Superspace 

Supersymmetry transformations are generated by operators Q which transform bosons into 
fermions, i.e. these generators have fermionic character p6|j27| ]. The operators Q and their Hermi- 
tian adjoints Q can be chosen to be Weyl spinorsQ, which obey anticommutation relations 

{Qa,Qp} = {Q„,Qp} = 0, (2.1) 

{Qa,Q a } = 2a a ^P fl , (2.2) 

[Q a ,P„] = [Q^P„] = 0, (2.3) 

where P^ is the energy-momentum operator. Together with the familiar commutation relations for 
the generators P^ and M^ u of the Poincare group, Q and Q form a closed algebra, the so-called 
super-Poincare algebra. 

A compact technique for working out representations of the supersymmetry algebra was pro- 
posed by Salam and Strathdee [28]. They introduced Grassmann variables 6 a (a = 1,2) and 9^ 



1 For our notations and conventions see app. |B|. 
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(a = 1,2), 



{e a ,9( , } = {o",o fi } = {6 a ,e f> } = o. (2.4) 



With these anticommuting parameters the supersymmetry algebra can be expressed in terms of 
commutators, 

[9Q,0Q] = [Q9,Q9] = , (2.5) 
[9Q,Q6] = WaKePp , (2.6) 
[9Q,P„) = [Q9,P I ,) = 0. (2.7) 

Now we are able to exponentiate the super-Poincare algebra into a group in such a way that the 
product of two group elements is again a group element. An element of this super-Poincare group 
is given by 

G(x,0,6) = j&P+OQ+OQ) . (2.8) 

These group elements generate transformations in the eight-dimensional superspace parametrized 
by the coordinates (x^, 9 a , 9^). In the following we use z = (x^,9 a ,9^) to denote a point in 
superspace. 

Left action of the group element G(a, rj, rj ) induces a motion in superspace, 

(x^,9,9) (a^ + a M + vqa^O - i9a^f},9 + r},0 + rj) . (2.9) 
This transformation is realized by the differential operator representation of the algebra 



P» = idp , (2.10) 
8_ 



Qa = ^+i^9 a d tl , (2.11) 



Q^ = — 3 — - i9 a (7-^d u , (2.12) 
89 



where differentiation with respect to spinor parameters 9 and 9 is defined by 



3 -y = % , = 8{ . (2.13) 

89 c 



89" " 



The usual rules for raising and lowering spinor indices (cf. app. |B|) therefore give an additional sign, 
if the index position in the differentiations is changed, 
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Correspondingly, right action of the group elements induces an anti-realization of the super- 
Poincare group generated by covariant derivatives D and D, 

^ = ~^%, (2-15) 
»a = 4 + (2.16) 



de a 



which obey the anticommutation rules 



{D Q ,D /3 } = {D^,D. /3 } = (2.17) 
{D a ,D.} = 2a a .»P,. (2.18) 

Furthermore, one can derive the following useful identities 

D a D (3 = _}_ £ a(3 D 2 ^ ^.19) 

D^D 13 = ^e^D 2 , (2.20) 

D a D 2 D a = D ii D 2 D a , (2.21) 
D 2 D 2 D 2 = —16 nT) 2 , (2.22) 
D 2 D 2 D 2 = -16 nD 2 , (2.23) 

where the squared covariant derivatives are given by 

D 2 = D«D a = -e"^^ + 2i(9v^y» + 9 2 n, (2.24) 

d 2 = d^d* = d _. 9 . +2i(ea' x -^)d lx + e 2 n . (2.25) 
2.2 Superfields 

Supersymmetric theories are most easily formulated in terms of superfields in superspace. In order 
to get a feeling for how to define a superfield, let us first consider an ordinary quantum field (j)(x) 
which depends only on the coordinates of Minkowski space. Translations of these coordinates 
are generated by the operator P^, and we can consider <p(x) to have been translated from x^ = 0, 

<f,( x ) = e lx - p 4>(0)e- ix - p . (2.26) 
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In complete analogy, a superfield F(x,9,9) can be denned as 27] 



F(x, 9,9)= G(x, 9, 8 )F(0, 0, OjG" 1 (x, 9,9) , (2.27) 



where G(x,9,9) is an element of the super- Poincare group given by eq. ( |2.8|) . This means that a 
superfield is defined as a Taylor expansion in 9 and 9 with coefficients which are themselves local 
fields in Minkowski space. Due to the Grassmann nature of 9 and 9, this expansion breaks off, and 
the most general superfield reads 

F{x,9,9) = f(x) +9(^{x) +9x(x) + 9 2 m(x) +~9~ 2 n(x) +0a»0v IJt {x) 

+9 2 9X(x) + 9 2 9il>{x) + 9 2 9 2 d(x) . (2.28) 



This superfield contains as Taylor coefficients four complex scalar fields /, m, n and d, one complex 
vector v^, two spinors <p an d ^ i n the (|,0) representation and two spinors x an d A in the (0, ^) 
representation of the Lorentz group, altogether 16 fermionic and 16 bosonic field components. 

Consequently, superfields form linear representations of the supersymmetry algebra which are, in 
general, highly reducible. Irreducible representations can be constructed by imposing constraints on 
the superfields. Like all covariant derivatives, D and D can be used to impose covariant conditions. 
Chiral superfields $ are characterized by the condition 

= . (2.29) 

This first order differential equation is most easily solved in terms of the variables = x^ — i9a^9 
and 9 since 

(x" - i9a^9) = and DJ& = . (2.30) 
Then an arbitrary function of y and 9 is a chiral superfield, 

<%, 9) = A{y) + V29i,{y) + 9 2 F{y) . (2.31) 

This is the most general solution to eq. (2.29), as may be seen by expressing the covariant derivatives 



in terms of y, 9 and 9, 
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By Taylor expansion in 9 and 9, we can write a chiral superfield as a function of the original 
superspace coordinates x^, 9 and 9, 



$(x,6,6) = A(x) -i9a< 1 9d fl A(x) -^9 2 9 2 nA(x) 

+V29tp(x) + ■^ r 9 2 d^{x)a il 9 + 9 2 F{x) . 
v2 

Conjugation gives an antichiral superfield $ which satisfies the constraint 

Av¥ = . 



(2.34) 



(2.35) 



It is a natural function of = x^ + i9a^9 and 9, and its power series expansio is obtained from 
eqs. ( 2.31 ) and (|2.34 ) by conjugation. 

Supersymmetry invariant actions can be constructed from chiral superfields and their products. 



It is clear from the expansion (2.31) that a product of chiral superfields is again a chiral superfield 
(cf. app. |B.3| ), whereas a product of a chiral and an antichiral superfield will satisfy neither eq. ( [2.29D 
nor eq. ( 2.35|) . However, not every component of these product superfields can be used to construct 
supersymmetric actions. To be able to formulate supersymmetric theories we have to isolate the 
components which are invariant under supersymmetry transformations, up to total derivatives. 



2.3 Superspace Invariants 

The general method by which a translation invariant action is derived from fields is to integrate 
a Lagrange density £(x) over d 4 x. The result is translationally invariant if surface terms vanish. 
Similarly, SUSY invariant actions can be constructed by integration over superspace, once we have 
defined an integral over the Grassmann variables 9 and 9. This Berezin integral is determined by 
imposing linearity and translation invariance, except for the normalization which is fixed by the 
definitions p9[] 

jd 2 99 2 = l and Jd 2 9 9 2 = l, (2.36) 
with all other integrals vanishing. The two-dimensional volume elements are defined by 

d 2 9 = ~d9 a d9 a , (2.37) 

d 2 9 = --d0^d0 a . (2.38) 
For integration over superspace we introduce the following integration measures 

d 6 s = d 4 x d 2 9 , d 6 s = d A x d 2 9 , (2.39) 
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d 8 z = d^xd l 9d l 9. (2.40) 

If we adopt the convention of dropping total divergences, i.e. surface integrals, the differential 
operators — \D 2 and — \P> 2 are equivalent to d 2 9 and d 2 6 under a space-time volume integral, 

J d 8 zF(x,9,9) = J d 6 s (-^) F {x,9,9) ( 2 - 41 ) 

= J f d 6 s (-^ 2 ) F(x,0,0) (2.42) 

= J d A x^-F(x,e,d) (2.43) 

= Jd 4 x?^F(x,9,e), (2.44) 
where F is an arbitrary function of x, 9 and 9, i.e. a superfield. It follows that 

J d 8 zD a F(x,9,9) = J d 8 zD^F(x,9,6) =0 . (2.45) 
Hence, we immediately get the following rules of integration by parts 

J d 8 z F 1 D a F 2 = ^J d 8 z {D a F 1 )F 2 , (2.46) 

J d 8 z F 1 D a F 2 = T J d 8 z {D a Fi) F 2 , (2.47) 

where the upper (lower) sign is valid if F\ is an even (odd) Grassmann function. Similarly, higher 
powers of covariant derivatives can be partially integrated by means of the following formulae 

Jd 8 zF 1 (D 2 F 2 ) = J d 8 z (D 2 F 1 ) F 2 , (2.48) 

J d 8 zF 1 (D 2 F 2 ) = J d 8 z (D 2 F^) F 2 , (2.49) 

J d 8 zF 1 (p 2 D 2 F 2 ^j = j d 8 z (D 2 D 2 F^ F 2 , (2.50) 

J d 8 zF 1 (D a D 2 D a F 2 ^ = =f J d 8 z (d 2 D a F^j D a F 2 , (2.51) 

= J d 8 z (D a D 2 D a F^j F 2 . (2.52) 
We may also define superspace delta distributions 

5 2 {9-9') = {6-9') 2 , 5 2 {9-W) = (9-W) 2 , (2.53) 

5 s (z-z') = 5\x - x')5 2 {9 - 9')5 2 {9 -W) . (2.54) 



2.3. SUPERSPACE INVARIANTS 



23 



Applying covariant derivatives to these delta functions yields 

D 2 5 2 (9 1 - 9 2 ) = -4exp[i(ei-e 2 )a^d ltfl ] , (2.55) 

D{ 2 P(0l-h) = -4exp[-^a^(^-^)9i iM ] , (2.56) 

TT 1 2 D 2 5 2 {9 1 - 9 2 )5 2 (9~ 1 -T 2 ) = l6exp[-i(e 1 a»e~ 1 + 8 2 aL%-2e 1 a»6 2 ~)d 1 ^] , (2.57) 

D 2 TT 1 2 5 2 {e 1 - 6 2 )5 2 {T 1 -T 2 ) = 16exp[i(6 1 ^% + e 2 a^e 2 :- 29 2 ^%)d l4l ] , (2.58) 

where D\ and D\ act on x%, 9\ and Q\. The argument of the derivatives can be changed through 
the following transfer rules 

= -D 2 *5 8 (z 1 -z 2 ) , (2.59) 

= -TT 2 ^8 8 (z 1 -z 2 ) , (2.60) 



DfS 8 


>i 


- z 2 




>i 


- z 2 


D 2 5 8 


>i 


- z 2 


TT^d 8 


>i 


- z 2 



(2.61) 

(2.62) 



and since covariant derivatives with different arguments anticommute we also have 



TT l 2 D 2 1 5 8 {z 1 -z 2 ) = TT 1 2 Dl5 8 {z 1 - z 2 ) (2.63) 
= D 2 TT 1 2 5 8 {z 1 -z 2 ) 
= D 2 TT 2 2 5 8 { Zl -z 2 ). 

Superspace integration can be used to construct invariant actions. Consider first the integral 
over a chiral superfield. Due to the constraint ( 2.2S| ), a chiral superfield is independent of 9, i.e. 



J d 2 9, and hence the full superspace integral gives zero. 

Since for chiral superfields the supersymmetry algebra can be realized as coordinate transfor- 
mations of the chiral subspace of superspace alone, which has coordinates and 9 a but not 9- , 
the J d 4 xd 2 9 integral, without the d 2 9, is already an invariant integral for chiral superfields. There- 
fore, the most general supersymmetric renormalizable Lagrange density involving only one chiral 
superfield reads 



C = I d 2 9d 2 9$$ + 



d 2 9 ( ^m$ 2 + ^A$ 3 + g§ ) + li.c. 



(2.64) 



We will omit the tadpole term g$> in the following, since it can always be eliminated by field 
redefinitions. 
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2.4 Superfield Propagator 



In close analogy to the usual perturbation theory one can develop a perturbation theory in super- 



space [ 30 , 25 1 . Our goal is to compute Green functions for superfields, 



(z\...,z r ;z r+ \...,z N ) = (0\T{&(z L )---Q(z r )$(z r+1 )~-$(z N )}\0) , (2.65) 



where z l = (x lM , 9 ia , 9 i -) denotes a point in configuration superspace. Let us start by computing 



the propagator for a chiral superfield, which is constructed from the free Lagrangian 



Ln 



ft.: <]><!> -- ( / (fs ^m$ 2 + h.c. 



(2.66) 



Since the operator — D D 2 /(16 □) projects on chiral fields, 

1 D 2 D 2 - 

-<!» = $ if D$> = , 



16 □ - - < 2 ' 67 > 

we can use eqs. ( 2.41[ )-( p.44 ), and rewrite the d 6 s integration in the mass term of the free Lagrangian 
into an integration over the whole superspace, 



1 



D 2 



-D 2 



drz < + -m $ — $ + $ $ 



□ 



□ 



(2.68) 



with the matrix 



/ I 



m 



M 



-—D £ 

4n 



V 



1 



1 \ 



--D 2 j 

4n / 



(2.69) 



where we have assumed a real mass m. 

To derive equations of motion we have to define a functional derivative in superspace, where 
we have to take into account the chirality constraint D $ = 0. This constraint is automatically 
respected by varying in the y basis, 

8 



<S>{y\9')=5\y-y')5 2 (9-9') 



(2.70) 



Going back to the variable x, the variation under superspace integration reads 



tf*(a?,0,0) 



d 8 z' $(x',9',0>)F(x',9',9 r 



\d 2 F(x,9,9[ 



(2.71) 
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This leads to the formal definition 

6 



-$( x ',e',e') = -\d 2 5 8 (z- z ') 
5$(x,e,e) K 1 4 v ' 



Variation of the free Lagrangian ( 2.68j ) then gives equations of motion, 

. 



1 / D 2 



M ( t 



D 



With eq. (2.67) this leads to 



m& - -D 2 
4 

— 1 9 

m$ D $ 

4 



0, 
. 



The propagator is defined as Green function of the operator 



1 / D 

D 2 



M 



i.e. the differential equation defining the two point function A reads 

MA = 



1 I D 2 



4 □ 

D~ I n 1D2 

\ TE J 

where the differential operator on the right-hand side implements the chirality constraint ( |2,29| ), 
Solving this equation, one gets the propagator for a chiral superfield 



5{z - z') 



(2.72) 
(2.73) 

(2.74) 
(2.75) 

(2.76) 



(2.77) 



A(z,z') 



1 



□ + m 1 



1 „ 9 ^2 



16 



D D z 



m 



— D D —D z 



\ 16 



8(z - z') 



(2.78) 



where — !/(□ + n, 2 ) is a symbolic notation for the Green function of the Klein-Gordon operator 



□ + 77T 



2.5 The Generating Functional 



The generating functional for free Green functions is given by the vacuum-to- vacuum amplitude in 
the presence of an external classical chiral source J coupled to a free chiral field 

/ 1D " n \ 
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(2.80) 



where Lq is the free Lagrangian from eq. ( 2.68j ), and N is a normalization factor which can be chosen 
such that Zq[0, 0] = 1. In non-super symmetric quantum field theories the role of this normalization 
factor is to take out disconnected vacuum bubbles, which would otherwise contribute to Green 
functions. Although vacuum diagrams vanish in super symmetric theories because of the non- 
renormalization theorems (cf. section [2.6) ), the normalization factor does not equal unity when the 
volume of the system tends to infinity 

The oscillatory path integral fl2.80| ) is not well defined, and has to be Wick rotated to Euclidean 
space to be evaluated unambiguously. The Green functions calculated in Euclidean space then 
yield Green functions in Minkowski space by analytic continuation. We will write all quantities in 
Minkowski space with the understanding that they can be justified in Euclidean space. 

Performing the functional integral with standard techniques, one gets 



Zq [J, J] = exp < 



- / d*zd 8 z' (J(z), J(z)) 



( 1^! 

4 □ 







\ 



V 



o 



4 □ / 



A(z,z') x 



1~ 







\ 



V 







ID' 

1~n ) 




(2.81) 



where A(z,z') is the chiral superfield propagator ( 2.78 ). With eq. ( 2.67 ) this can be brought to a 
familiar form, 
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(2.82) 



where A GRS (z,z') is the superfield propagator of Grisaru, Rocek and Siegel |32], 
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Since J is a chiral source, its functional derivative is defined like in eq. fl2.72| ). The functional 
derivative of Zq then reads 



/ 1 i \ 

i 6J(z) 

1 5 



Zq 



dV A(x,a/) 



4n 



/ D 2 J{z 

\ i SJ(z) I v 
With eqs. ( p. 77 ) and fl2.67| ) this yields a functional equation for Zq, 
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\ iSJ(z) J 

This equation can easily be generalized to the interacting case. For the <3? 3 theory coupled to an 



external source the equations of motion read 



(2.8 



By comparison with the functional equation for Zq ( 2.85 ) we can write down the defining equation 
for the full generating functional Z[J, J], 
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Using the interaction Lagrangian £i NT 

Ant = J d 2 9^ 3 + J d 2 9h$ 3 , 

the operator on the right-hand side of eq. ( |2.87| ) can be rewritten as 
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Hence, eq. ( p. 87 ) yields 
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By comparison with the functional equation for Zq ( |2.85| ), the generating functional Z can be 
related to the free-field generating functional 



Z[J, J] = e 



18 18 
- INT V i SJ' i 8 J 



Z [J, J] . 



(2.91) 



This relation, familiar from ordinary quantum field theory, is the starting point of perturbation 
theory in superspace. 



2.6 Feynman Rules in Superspace 



iV-point Green functions are obtained from the generating functional by functional derivation, 
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where we have used eq. ( |2.91 ) in the last step. The factors 



U, X A-TNT I • r T ) • ^ 

I 6J I SJ 



(2.92) 



(2.93) 



(2.94) 



generate vertices at the superspace points z' J , and the functional derivatives jj in £ INT , when acting 
on Zq, generate propagators connecting different vertices. The operators jjj-rj riot in £i NT generate 
propagators on external lines, which have to be amputated and replaced by superfields $>(z l ) in 
order to get an effective contribution to the Lagrangian. This leads to the following Feynman rules 
in configuration superspace 
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1. Provide external lines with factors 3>(z) or $(z). 



2. Due to the chirality of the functional derivative (cf. eq. ( 2.72; )), one has to include one (two) 



factors —\D 2 acting on internal propagators to each $ 3 vertex with two (three) internal lines. 
The same applies to <J> 3 vertices and factors — jD 2 . 

3. Include a coupling constant |A and a superspace integration J d 8 z for each vertex. 

4. Use the GRS-propagators iA GRS ( 2.83| ) for internal lines. 



5. Include the usual symmetry factors. 

In the next chapter we will illustrate these Feynman rules by performing several sample calculations. 
Let us first investigate the general structure of diagrams calculated with these rules. 

The d 2 9 l d 2 9 % integrals at each vertex can be done, leaving us with one overall d 2 9 d 2 9 integral 
and the usual Minkowski space integrals over d x l for each diagram. To see how this comes about let 
us follow the ^-integrations around an arbitrary closed loopQ. It consists of propagators, including 
factors S 8 (zi — Zj+i) and covariant derivatives acting on them, external superfield factors, and 
d 2 9 l d 2 9 % integrals. Higher powers of covariant derivatives can be reduced by using the identities 

0-(PD. 

Consider now the propagator from one vertex z 1 to another one z 3 , and integrate by parts 

using eqs. ( |2.46| )-( |2.52j ) to remove all the covariant derivatives from its 5-function. The original 
contribution becomes a sum of terms. If there are other propagators connecting z % and z 3 we can 
use the relations 

SijSij = , (2.95) 

5 ij D a S ij = , (2.96) 

5 ij D 2 5 lj = , (2.97) 

5 ij D a D a 5 ij = , (2.98) 

5 lj D a D 2 5 ij = , (2.99) 

5 tj D 2 D 2 5 ij = 5 lj D 2 D 2 5 ij = 6 ij D a D 2 D a 5 ij = 5^ D^D 2 D a <% = 164,- <5 4 (x; - Xj ) , (2.100) 

5ijD a D 2 D@5ij = Se^Sij 5 4 ( Xi - Xj ) , (2.101) 

8 ij D a D 2 D l3 S ij = 8e?hij 5 4 ( Xi - Xj ) , (2.102) 



2 We will assume that loop divergences have been properly regularized. 
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where Sy = 5 8 (zi — Zj). Hence, the terms generated by the partial integration vanish, unless each of 
the other ^-functions has exactly two D's and two -D's acting on it. Now the free <5-function can be 
used to perform the d 2 0i d 2 9 3 integral and shrink all the propagators between z % and zi to a point 
in #-space. This procedure can be repeated, until we have removed all 5-functions and performed 
all ^-integrals except the original one at z l . We are left with a sum of terms, all with a single d 2 9 d 2 9 
integral, various d'Alembert operators from eqs. ( 2.22 ) and ( 2.23| ), as well as covariant derivatives 
acting on the external superfields. 

Hence, we have ended up with a d 2 9 d 2 9 integral, even though in the original Lagrangian we 
may have had chiral d 2 8 integrals. This is the perturbative no-renormalization theorem for chiral 
superfields [p3|| : radiative corrections do not induce renormalizations of F-terms, i.e. purely chiral 
mass or interaction terms. Furthermore, all vacuum diagrams vanish, since the d 2 6d 2 9 integral 
without any external super field vanishes. 



Chapter 3 

Supersymmetric Leptogenesis 



In this chapter, which is based on ref. 17], we present the supersymmetric generalization of the 
leptogenesis scenario suggested by Fukugita and Yanagida [12]. After having introduced the super- 
potential, we compute all the relevant decay widths, CP asymmetries and scattering cross sections. 
In order to check the results we have performed two independent calculations. First by using 
the component field Feynman rules from appendix Q and then by using the superfield techniques 
introduced in the last chapter. 



3.1 The Superpotential 

In supersymmetric unification scenarios based on SO(10), the effective theory below the (B — 
L) breaking scale is the MSSM supplemented by right-handed Majorana neutrinos. Neglecting 
soft breaking terms, the masses and Yukawa couplings relevant for leptogenesis are given by the 
superpotential 

W = ^N C MN C + ^H ie H 2 + HxeQ\ d D c + H x eL\E c + H 2 eQX u U c + H 2 eLX u N c , (3.1) 

where we have chosen a basis in which the Majorana mass matrix M and the Yukawa coupling 
matrices A^ and A/ for the down-type quarks and the charged leptons are diagonal with real and 
positive eigenvalues. The corresponding Lagrange density reads 

£ = Jd 2 9W + J d 2 6W . (3.2) 

The chiral superfields in the superpotential ( |3.1| ) are most conveniently parametrized in the 
y-basis (cf . section |2.2| ) . 

Hi(y,0) = H t (y) + V2eHi(y) + 9 2 F H .(y) , (3.3) 
31 
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Q(y,9 
L(y,6 
U c (y,9 
D c (y,9 
E c (y,9 
N c (y,6 



q(y) + V2eq^y)+e 2 F Q (y) , 
T(y) + V2ei L (y) + e 2 F l (y) , 
U^y) + V2eu R c (y) + 9 2 F u o(y) , 
rr(y) + V2ed B c (y)+9 2 F D c(y) , 
Wc(y) + V2 0e R c (y)+0 2 F EC (y) , 
N-(y)+V26v^{y) + e 2 F N c{y) . 



(3.4) 
(3.5) 
(3.6) 
(3.7) 
(3.8) 
(3.9) 



Q and L stand for the left-handed quark and lepton doublets, U c , D c , E c and N c are the right- 
handed singlet fields, and Hi denotes the two Higgs-doublets, 



Hi 



-Hi 



and 



Ho 



HI 
Hi 



(3.10) 



Besides the usual bispinors for quarks and charged leptons we can introduce Majorana-spinors 
for the right- and left-handed neutrinos, 



and v 



(3.H) 



In the symmetric phase of the MSSM no mixing occurs between the fermionic partners of gauge 
and Higgs bosons. Therefore, we have two Dirac higgsinos 



h° 



I H® 

1 a 



and h 



H{ 



Hi 



(3.12) 



\H° 

which again form an isospin doublet, 



h° 
-hr 

The auxiliary component fields of the chiral superfields ( |3.3| ) 
grange density |3l], 

'dc_\ ] 

dFj ■ 

The superpotential (|3.1| ) then yields the following contributions to the auxiliary fields 



F n 



nHl+[q\ d D c ) +{l\ l E 



(3.13) 

are obtained from the La- 
(3.14) 

(3.15) 
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H 2 



HH\ + (qX u U c ) + (l \ v N C J e , (3.16) 

(3.17) 
(3.18) 
(3.19) 
(3.20) 
(3.21) 
(3.22) 

The vacuum expectation values of the neutral Higgs fields generate Dirac masses for the down- 
type quarks and the charged leptons 
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'Hluf {\{) ji + H\Dj 
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F Ef 
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MiNf + eH\ 





vi = (H^) / m d = X d v 1 and mi = Xi vt , 

and for the up-type quarks and the neutrinos 



V-2 



m u = X u v 2 and m D = X u v 2 ■ 



(3.23) 



(3.24) 



The Major ana masses M for the right-handed neutrinos, which have to be much larger than the 
Dirac masses m D , offer a natural explanation for the smallness of the light neutrino masses via the 



see-saw mechanism [11]. 

To generate a non-vanishing baryon asymmetry, one needs a hierarchy in the Majorana mass 
matrix M. Then the scale at which the asymmetry is generated is given by the mass Mi of the 
lightest right-handed neutrino. Hence, it is convenient to write all the masses and energies in units 
of Mi, 

( M i\ 2 s , Mi . 

0"> = tt~ i x = — n and z = — — , (3.25) 
3 \Mi) T v ; 

where Mj are the masses of the heavier right-handed neutrinos, s is the squared centre of mass 

energy of a scattering process and T is the temperature. 



3.2 The Decay Channels of Heavy Neutrinos 

Since at these energy scales supersymmetry breaking terms can be safely neglected, we are working 
in a theory with exact supersymmetry, i.e. we can use the superfield techniques exposed in the 
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(a) 




Li 



Ho 



+ m 




N c N c + JV? 




Figure 3.1: Decay modes of the right-handed neutrino superfield. 



last chapter to compute decay and scattering amplitudes. To leading order, the decay modes of 
the right-handed neutrino and its scalar partner with positive lepton number in the final state 



are all contained in the superfield diagram fig. 3.1a, corresponding to the superpotential term 
—N c \t Lei/2- Choosing the usual decomposition of the 5-matrix 



S=l + iT , 

we get the following tree level contribution to the T-matrix 



(3.26) 



d 4 x / d 2 9N c \lLeH 2 



(3.27) 



The decay amplitudes we are interested in are then given by the matrix elements of iT^ . Let 
us first concentrate on two-particle final states. According to the component field decompositions 
( |3.3[ )-( |3T9| ) of the chiral superfields, the right-handed Majorana neutrinos Nj can decay into a lepton 
and a Higgs-boson or into a slepton and a higgsino, while their scalar partners N? can decay into 



a lepton and a higgsino or into a slepton and a Higgs boson (cf. fig. 3.2). The decay widths at tree 



level read [18| 
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(3.28) 



According to eq. ( D.15| ), the reaction densities for these decays are then given by 

(2) _ M? (mUfl)jj a j v / «7 1 



7iv, 



2 7 i 



-Kx{zy/cu) 



(3.29) 



(3.30) 
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Figure 3.2: Decay modes of the right-handed Majorana neutrinos and their scalar partners. 
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All these decay modes are CP violating, the dominant contribution to CP violation coming 



about through interference between the tree level and the one-loop diagrams shown in fig. 3.2. 



These one-loop diagrams can be summarized by the one-loop superfield diagrams in figs. |3.1| b and 
|3.l| c. Note that the non-renormalization theorems do not apply to these diagrams, since they are 
contributions to the D-term NjLieH 2 . Using the superspace Feynman rules discussed in section 



the diagram in fig. |3.l| b yields the following contribution to the T-matrix 
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After partial integration of the D a {z 3 ) derivatives, the only non-vanishing contribution reads 
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(3.33) 



where we have used eqs. (2.61), ( 2.62| ) and (2.23) in the last step. Now we can perform the d 2 9 3 d 2 9 3 
integration, 
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(3.34) 



Similarly, we can partially integrate the D 2 (z 2 ) derivatives, and perform the d 2 9 2 d 2 9 2 integrations 
after having used eq. (2.100) to remove all covariant derivatives. Then iT^ reads 



iT 
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(3.35) 



3.2. THE DECAY CHANNELS OF HEAVY NEUTRINOS 



37 



By Fourier transforming the loop propagators, we get 



iT, 
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At 
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x / d 4 xi d 4 x 2 d 4 x 3 
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x J d 2 9d 2 9N^(x 1 ,9,9)L i (x 2 ,9,9)eH 2 (x 3 ,9,9) . (3.36) 

Analogously, we can compute the contribution of the diagram in fig. |3.1| c to the T-matrix, 
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Again using the Fourier representation for the propagators, we find 
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(3.38) 



where Bq(— k,0, 0) is the massless two-point scalar integral defined in eq. ( |A.5| ). One-loop correc- 
tions to the decay amplitudes are given by the matrix elements of iTfp + iT^ . 

Interference between these one- loop diagrams and the tree- level amplitudes gives rise to CP 
asymmetries in the different decay channels of Nj and iV?, which can all be expressed by the same 
CP violation parameter e,-, 
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(3.39) 



x — 1 



for x 1 . 



Here n is the flavour index of the heavy (s)neutrino in the loop. This result agrees with the one in 



ref. [18 1 and is of the same order as the CP asymmetry in ref. [16|. 
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Figure 3.3: Contributions of the scalar potential to the decay width and the interactions of a scalar 
neutrino. 



With £j we can parametrize the reaction densities for the decays and inverse decays in the 
following way 
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(3.40) 
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Additionally, the scalar potential contains quartic scalar couplings, which enable the decay of 
Nj into three particles via the diagram shown in fig. |3.3| a. This is just the contribution of the 
auxiliary field F H given in eq. Q3.16 ) to the superfield decay amplitude in fig. [Q| a. The partial 
width for this decay is given by 
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and the corresponding reaction density reads 
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(3.44) 



(3.45) 



Since the Yukawa coupling of the top quark and its scalar partner is large, a u can be of order one. 

(3) (2) 

But even then j~ is much smaller than 7 ~ . Hence, the three particle decays give only a small 



3.2. THE DECAY CHANNELS OF HEAVY NEUTRINOS 



39 




N C N C 

2l> <Z2 




Li 



Ho 



Figure 3.4: Lepton number violating scatterings mediated by a neutrino superfield. 



correction, which we have taken into account for completeness. However, we have neglected the 



CP asymmetry in this decay which comes about through the one-loop diagrams in fig. 3.1 



Furthermore, we have neglected the leptonic auxiliary field F L given in eq. ( 3.18 ), since its 
contribution to the sneutrino decay width will be of order (Aj.A 2 Aj,)jj, i.e. much smaller than the 
other partial decay widths, at least for the lightest sneutrino (j = 1). 

The dimensionless squared total decay widths of Nj and iV? are then finally given by 
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(3.46) 



(3.47) 



The vertex in fig. [Jl|a also gives 2^2 scattering processes involving one scalar neutrino, like 
N? + I —> q + U c (cf. fig. |3.3| b). The reduced cross section for this process reads 



<5"22, (x) = 3a u 



(m D m D )jj x — aj 



(3.48) 



For the processes iVJ + 



V + U c and Nf + U 



rt 



I + q, the corresponding back reactions 
and the CP conjugated processes we find the same result. The corresponding reaction density can 
then be calculated according to eq. ( p,16| ). One finds 

3a u M? (m) D m D )jj 



722, (z) 



16vr 4 



-3-Ki(zv^7) 



3a M 

47T QjZ 2 



In, (z) • 



(3.49) 



f2) > 

Hence, 722^ will be much larger than 7^. and 7^ for small ajZ 2 , i.e. for high temperatures T 3> Mj. 



Together with similar scatterings, which we are going to discuss in section |3,4| , these processes will 
therefore be very effective in bringing the heavy (s)neutrinos into thermal equilibrium at high 
temperatures where decays and inverse decays are suppressed by a time dilatation factor. 
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Figure 3.5: L violating processes mediated by a virtual Majorana neutrino or its scalar partner. 

3.3 Lepton Number Violating Scatterings Mediated by Right- 
Handed Neutrinos 



Using the tree level vertices from figs. |3.2| and |3.3| as building blocks we can construct lepton number 
violating scatterings mediated by a heavy (s)neutrino. Although of higher order than the tree level 
decays, these diagrams have to be taken into consideration to avoid the generation of an asymmetry 



in thermal equilibrium, which is forbidden by CPT invar iance [35,36]. 

All these processes can be summarized in the configuration space superfield diagram shown in 
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Figure 3.6: Diagrams contributing to the lepton number violating scatterings via heavy sneutrino 
exchange. 



fig. |3.4| . It yields the following contribution to the T-matrix 

^4^, e ik(x 1 -x 2 ) 

k 



x [ d 2 9d 2 eL J (x 1 ,e,e)eH 2 (x 1 ,e,9)^^L i (x 2 ,e,e)eH 2 {x 2 ,d,e) . 
J 4 n 2 



The superfield product can be evaluated with eq. (B.56). The matrix elements of iT AL then cor- 
respond to the component field processes that we are going to discuss in the following. In this 
section we will only mention the different processes which have to be considered. The correspond- 
ing reduced cross sections can be found in appendix [E] and the reaction densities are discussed in 
appendix [F]. 



By combining two of the decay vertices (cf. fig. ^2 and fig. |3.3| a) one gets the processes that we 
have shown in fig. |3.5| and the corresponding CP conjugated processes. We will use the following 
abbreviations for the reaction densities 

if = l(j+ H 2 <-» V + + q) , 7 ^ 6) = 1 (l + H 2 ^T+h 

if =^(l + h~T+ q^ + U cj 
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Figure 3.7: L violating processes mediated by a right-handed (s)neutrino in the t-channel. 



The contributions from on-shell (s)neutrinos contained in these reactions have already been taken 
into account as inverse decay followed by a decay. Hence, one has to subtract the contributions 
from real intermediate states to avoid a double counting of reactions [j36||. 



From the scattering vertex in fig. 3.2 b and the decay vertices we can construct the following 
processes 
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-f[l j +q^l + h + U cj ) =-f(l j + U C <->l + h + q1 
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In fig. 3.6 we have shown one typical diagram for each of these reaction densities. Again, these 



diagrams have on-shell contributions which have to be subtracted, since they can be described as 
decay of a sneutrino which has been produced in a scattering process. 

Up to now we have only considered processes with a neutrino or its scalar partner in the s- 



channel. In fig. |3.7| we have shown a selection of diagrams without on-shell contributions. The 
corresponding reaction densities will be denoted by 
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T+ I *-> h + H\) , 

H 2 + gt «_ f + P + fjc) = 1 {h 2 + uc ] ^ J 1 + V + 
T+l++h + H 2 ) , 7i 18) =l(j+Hl^h + l 

l+V ^h + p + lT c1 ) =j(l + q^T+u' c1 + h 
l + V c ^T+q ] +h) =j(T j + h^~l + q j + JJ J 



+ h^l + l + U c =j[U c + h^l + l + q J ) . 



At first sight one may think that these diagrams could be neglected, since they are suppressed 
at intermediate temperatures, i.e. intermediate energies x ~ aj. However, they give an important 
contribution to the effective lepton number violating interactions at low energies and therefore have 
to be taken into consideration. 



3.4 Interactions with a Top or a Stop 

The Yukawa coupling of the top quark is large. Thus we have to take into account lepton number 
violating interactions of a right-handed neutrino with a top quark or its scalar partner. In addition 




Figure 3.8: Neutrino- (s)top scattering in configuration superspace. i, j, q andp are flavour indices. 
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Figure 3.9: Neutrino- (s)top scattering. 



to the processes already considered in section |T2| (cf. fig. pT^b), we have the superfield diagram 



shown in fig. 3.8, which gives the following contribution to the T-matrix 



iT, 



-i(K)ij(\i) pq / d 4 xid A x 2 



^4^, e ik(x 1 -x 2 ) 



(2vr) 4 k 2 + is 



x / d 2 ed 2 6N^x 1 ,e,9)(L i (x 1 ,6,e)Q q {x2,e i 6))U^x 2 ,8,\ 



(3.51) 



In component fields, we have the following processes with a Majorana neutrino Nj as external line 



(cf. fig. po 
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-yfo + U$ + q \ = j^Nj + q^V +u\ , 

j^Nj + l^q + ll] , 

"/(Nj + u++l + q) = 7 (Nj + q^l + u) . 



(4) 

At this order of perturbation theory these processes are CP invariant. Hence, we have the same 
reaction densities for the CP conjugated processes. 

For the scalar neutrinos we have similarly (cf. fig. 3. 10| ) 
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Figure 3.10: Sneutrino-(s)top scattering. 



l\Nf + l ^q + u 

l[N!f + q^T+u\ = jlN? + u<r+l + q 



(9) 
% 

The quartic scalar couplings from the scalar potential give additional 2 
processes, which can be neglected since they are phase space suppressed. 



3, 3 -> 3 and 2 4 



3.5 Neutrino Pair Creation and Annihilation 



The Yukawa couplings of the right-handed neutrinos also allow lepton number conserving processes 



like neutrino pair creation and annihilation. The two superfield diagrams in fig. |3.11| yield the 
following contributions to the T-matrix 
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Figure 3.11: Pair annihilation of singlet neutrino superfields in configuration superspace. 
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Figure 3.12: Neutrino pair annihilation. 
Decomposing the superfield product into component fields with the help of eq. QB.57| ), we get the 



processes depicted in fig. 3.12 for the neutrinos 



lN i N j 



jlNi + Nj^l + l 



J 2 ) 
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^L=l(Ni + Nj^H 2 + Hl 



For the scalar neutrinos we have similar diagrams and additional contributions from quartic 
scalar couplings (cf. fig. 3.13j ). We have the following reaction densities 



7- 



N9N9 
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It is interesting to note that the contributions to 7~~ and r y~~ from the scalar potential are 

& 'N C N C 'N9N c ^ 

%■ 3 * 3 



not contained in the superfield diagrams in fig. p,ll| . They originate from the contribution of the 
auxiliary field F H2 (cf. eq. |3.16 ) to the decay diagram in fig. jO|a, 



Finally, there are neutrino-sneutrino scattering processes (cf. fig. 3.14 ), 



%~ c =l[Nf + N 3 ~l + l ) , 7 W r r , N c + N .^ h + H2 

3 i \ / 3 i \ 

Such diagrams also give neutrino-sneutrino transitions like Nf + I «-> Nj + I. These processes 
transform neutrinos into sneutrinos and leptons into sleptons, i.e. they tend to balance out the 
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Figure 3.13: Sneutrino pair annihilation. 



number densities of the fermions and their super symmetric partners, but they cannot wash out any 
generated asymmetry. As we will see in the next chapter, the number densities of the neutrinos 
and the scalar neutrinos are already equal without taking into account these interactions, while 
the equality of the number densities of leptons and sleptons is ensured by MSSM-processes, which 
we are going to discuss in the next section. Finally, the dominant contributions to these neutrino- 
sneutrino transitions come from inverse decays, decays and scatterings off a (s)top which we have 
already considered. Hence, we can neglect these additional processes. 
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Figure 3.14: Neutrino- sneutrino scattering. 
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Figure 3.15: Example of a Lf and L s violating MSSM process. 

3.6 MSSM Processes 

In the MSSM the fermionic lepton number Lf and the lepton number stored in the scalar leptons 
L s are not separately conserved. There are processes transforming leptons into scalar leptons and 
vice versa. As an example we have considered the process e + e «-> e + e (cf . fig. 3.15| ) . For large 



temperatures, i.e. s 3> m~, the reduced cross section for this process is given by ||37| 

0"MSSM ~ 47TQ! 

This translates into the following reaction density 

7mssm 
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where we have introduced the dimensionless squared photino mass 

Mi 



(3.55) 



(3.56) 



These processes are in thermal equilibrium if the reaction rates are larger than the Hubble parameter 
H. This condition gives a very weak upper bound on the photino mass, 



™7 * 2.5 x 10 9 GeV { ^ ^ 



exp 



1 



412 V 10 10 GeV 



(3.57) 



In the calculations we assume m~ = 100 GeV. 



Chapter 4 



Numerical Results 



Now that we have identified all the relevant processes we can write down the network of Boltzmann 
equations which governs the time evolution of the neutrino and sneutrino number densities and of 
the lepton asymmetry 1 . In this chapter we work out the parameter dependence of the generated 
baryon asymmetry by solving the Boltzmann equations, and we discuss the role of the different 
scattering and decay processes |l7| . 



4.1 The Boltzmann Equations 



The evolution of the neutrino number Y N . as a function of the inverse dimensionless temperature 
z = Mi/T is given by 
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For the scalar neutrinos and their antiparticles it is convenient to use the sum and the difference 
of the particle numbers per comoving volume element as independent variables, 



(4.2) 



The Boltzmann equations for these quantities read 
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dz sH{M\) I I Y: 
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1 See app. ^ for a short review of kinetic theory in an expanding universe. 
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Furthermore, we have to discern the lepton asymmetry stored in the standard model particles Y Lf 
and the asymmetry Y Ls in the scalar leptons. Their evolution is governed by 
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+ ^ (7^ + 7 * L ) + ^ (7^ L " 7^ L ) + " p5U Tmssm (4.6) 




n, ^ n. 



y-eq ycq yeq yeq I 7 N N c I ) 

where we have introduced the following abbreviations for the lepton number violating scatterings 
mediated by a heavy (s)neutrino 

7T = 27«+7i 3) +7i 4) +7l 6) +7i 7) +27i r 12) + 7 r , (4.7) 

1% L = (4.8) 

7c Ai = 3 7 r +7 r+7i 18) +6 7 i 19) , (4.9) 

7 - = 4 7 ? + 27? + 8 7 ^ + 3 7 f + 4 7 i 15 ) + 2 7 f ) + 7 <v 17) + 7 £ 8 > + 6 7 ^ . (4.10) 

The numerical factors in front of the reaction densities arise due to the change in quantum numbers 
in the corresponding scattering, e.g. processes transforming leptons into sleptons appear with a 
relative minus sign in the Boltzmann equations for Y Lf and Y Lg . Furthermore, any reaction density 
is multiplied by the number of different processes (cf. chapter ||) contributing independently to the 
Boltzmann equations. 

This set of Boltzmann equations is valid for the most general case with arbitrary masses of 
the right-handed neutrinos. However, if the heavy neutrinos are mass degenerate, it is always 
possible to find a basis where the mass matrix M and the Yukawa matrix are diagonal, i.e. 
no asymmetry is generated. Therefore, one has to assume a mass hierarchy for the right-handed 
neutrinos, which in turn implies that the lepton number violating processes induced by the lightest 
right-handed neutrino are in thermal equilibrium as long as the temperature is higher than the mass 
of this neutrino. Hence, the lepton asymmetries generated in the decays of the heavier right-handed 
neutrinos are washed out and the asymmetry that we observe today must have been generated by 
the lightest right-handed neutrino. We will assume that the first generation neutrino N\ is the 
lightest. 

Hence, we will always neglect the heavier right-handed neutrinos as free particles. However, 
they have to be taken into account as intermediate states, since they give a substantial contribution 
to the effective lepton number violating processes at low energies. 
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The fermionic part Y Lf of the generated lepton asymmetry will be transformed into a (B — L) 
asymmetry by the action of sphalerons. But since MSSM processes like the one in section |3.6j 
enforce the relation 

Y Lf = Y La , (4.11) 
the total lepton asymmetry Yl = Y L + Y Ls will be proportional to the baryon asymmetry f38fl , 

Y B = -[ — — Yl , (4.12) 

where Nf is the number of quark-lepton families, and Nh the number of Higgs doublets. In our 
model with Nj = 3 and Nh = 2 we have 

Yb = ~^Y L . (4.13) 

From the observed baryon asymmetry 

Y B = (0.6 - 1) • 1(T 10 , (4.14) 



and eq. (4.11) we can infer the asymmetries that we have to generate, 

Y Lf = Y Ls = -(0.9 - 1.4) • 10- 10 . (4.15) 

The additional anomalous global symmetries in supersymmetric theories at high temperatures 
have no influence on these considerations, since they are broken well before the electroweak phase 



transition 3£]. 



4.2 The Generated Lepton Asymmetry 



Typical solutions of the Boltzmann equations are shown in fig. |4.1| , where we have assumed a 
neutrino mass M\ = 10 10 GeV, and a mass hierarchy of the form 

a 2 = 10 3 , {m ] D m D ) 2 2 = a 2 (m£,m )n , (4.16) 

a 3 = 10 6 , (mbm D )33 = a 3 (m[,m D )n . (4.17) 

Furthermore, we have assumed a CP asymmetry e\ = — 10 -6 . The only difference between both 
figures lies in the choice of {m) D m D )ii, 

„ (mU.jii / 10 ~ 4 eV for fi S- 41a ' 

m\ := — = < (4.18) 
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Figure 4.1: Typical solutions of the Boltzmann equations. The dashed line represents the equilibrium 
distribution for the neutrinos N\ and the solid lines show the solutions for the (s)neutrino number 
and the absolute values of the asymmetries in Lf and L s , while the dotted line represents the 
absolute value of the scalar neutrino asymmetry Yj_ . The lines for Y N and Y\ + and for the two 
asymmetries Y L and Y Ls cannot be distinguished, since they are lying one upon another. The 



hatched area shows the measured value (4-15). 



Finally, as starting condition we have assumed that all the number densities vanish at high tem- 
peratures T 3> Mi, including the neutrino numbers Y N and Y\ + . As one can see, the Yukawa 
interactions are strong enough to create a substantial number of neutrinos and scalar neutrinos in 



fig. 4.1a, even if Y N and Y\+ do not reach their equilibrium values as long as z < 1. However, the 
generated asymmetries 

Y Lf = Y Ls = -4 • 10- 10 (4.19) 



are of the requested magnitude. On the other hand, in fig. 4.1b the Yukawa interactions are much 
stronger, i.e. the neutrinos are driven into equilibrium rapidly at high temperatures. However, 
the large Yukawa couplings also increase the reaction rates for lepton number violating processes 
which can wash out a generated asymmetry, i.e. the final asymmetries are much smaller than in 
the previous case, 



Y T , 



Y, 



-6-10 



-12 



(4.20) 
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In both cases a small scalar neutrino asymmetry Yj_ is temporarily generated. However, Y%— is 
very small and has virtually no influence on the generated lepton asymmetries. 

Usually it is assumed that one has a thermal population of right-handed neutrinos at high 
temperatures which decay at very low temperatures T <C M\ where one can neglect lepton number 
violating scatterings. Then the generated lepton asymmetry is proportional to the CP asymmetry 
and the number of decaying neutrinos and sneutrinos @ , 

Y L « £1 [ Y%(T » Mi) + YlliT » M x ) ] ~ ^ • (4.21) 

With ei = — 10~ 6 this gives 

Y Lf = Y Ls « -2 • 1(T 9 . (4.22) 

By comparison with eqs. ( [4.19 ) and (|4.20D one sees that by assuming a thermal population of heavy 
neutrinos at high temperatures and neglecting the lepton number violating scatterings, one largely 
overestimates the generated lepton asymmetries. 

A characteristic feature of the non-super symmetric version of this baryogenesis mechanism is 
that the generated asymmetry does not depend on the neutrino mass M\ and {rrv D m D )ii separately 
but only on the ratio rh\ pi]. To check if this is also the case in the supersymmetric scenario we 



have varied m\ while keeping all the other parameters fixed. In fig. 4.2 we have plotted the 
total lepton asymmetry Yi = Y L + Y Ls as a function of rh\ for the right-handed neutrino masses 
Mi = 10 12 GeV, 10 10 GeV and 10 8 GeV and the CP asymmetry e x = -1(T 6 . 

The main difference between the supersymmetric and non-supersymmetric scenarios concerns 
the necessary production of the neutrinos at high temperatures. In the non-supersymmetric scenario 
the Yukawa interactions are too weak to account for this, i.e. additional interactions of the right- 
handed neutrinos have to be introduced. This is no longer the case here. The supersymmetric 
Yukawa interactions are much more important, and can produce a thermal population of right- 
handed neutrinos, i.e. the same vertices which are responsible for the generation of the asymmetry 
can also bring the neutrinos into thermal equilibrium at high temperatures. However, these lepton 
number violating processes will also erase a part of the generated asymmetry, hereby giving rise to 
the rh\ dependence of the generated asymmetry which we shall discuss in detail. 

First one sees that in the whole parameter range the generated asymmetry is much smaller 
than the naively expected value 4- 10 -9 . For low rh\ the reason being that the Yukawa interactions 



are too weak to bring the neutrinos into equilibrium at high temperatures, like in fig. [4,l| a. For 
high rh\ on the other hand, the lepton number violating scatterings wash out a large part of the 



generated asymmetry at temperatures T < M\, like in fig. 4.1b. Hence, the requested asymmetry 
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Figure 4.2: Generated {B — L) asymmetry as a function of mi for Mi = 10 8 GeV (dotted line), 
Mi = IO 10 GeV (solid line) and Mi = 10 12 GeV (dashed line). The shaded area shows the measured 
value for the asymmetry. 

can only be generated if fhi is larger than ~ 10 -5 eV and smaller than ~ 5 • 10~ 3 eV, depending 
on the heavy neutrino mass M%. 



The asymmetry in fig. 4.2 depends almost only on mi for small mi ^10 eV, since in this 



region of parameter space the asymmetry depends mostly on the number of neutrinos generated 
at high temperatures, i.e. on the strength of the processes in which a right-handed neutrino can 
be generated or annihilated. The dominant reactions are decays, inverse decays and scatterings 
with a (s)top, which all give contributions proportional to fhi to the Boltzmann equations at high 
temperatures, 



sH(Mi] 
— z 



7^ oc mi , 
722! oc fhi , 



sH(M_ 



Z (2) 

7~c oc mi , 



sH(Mi 



(3) 

lid oc mi 



7*i (r> Mi) oc fhi 



(4.23) 



sH(Mi) ,ZZ1 1 ' sH(Mi) 

For large fhi ^ 10 -4 eV on the other hand, the neutrinos reach thermal equilibrium at high 
temperatures, i.e. the generated asymmetry depends mostly on the influence of the lepton number 
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violating scatterings at temperatures T£,M\. In contrast to the relations of eq. ( [4.23| ) the lepton 
number violating processes mediated by a heavy neutrino behave like 

Tf^ocMx^mf, i = A,...,D (4.24) 



sH(Mi) 

at low temperatures. Hence, one expects that the generated asymmetry becomes smaller for growing 



neutrino mass Mi and this is exactly what one observes in fig. 4.2. 

Eq. (4.24) can also explain the small dependence of the asymmetry on the heavy neutrino mass 
Mi for mi ^> 1CP 4 eV. The inverse decay processes which take part in producing the neutrinos at 
high temperatures are CP violating, i.e. they generate a lepton asymmetry at high temperatures. 
Due to the interplay of inverse decay processes and lepton number violating 2 — > 2 scatterings 
this asymmetry has a different sign compared to the one generated in neutrino decays at low 
temperatures, i.e. the asymmetries will partially cancel each other, as one can see in the change 



of sign of the asymmetry in fig. 4.1a. This cancellation can only be avoided if the asymmetry 
generated at high temperatures is washed out before the neutrinos decay. At high temperatures 
the lepton number violating scatterings behave like 

7f L oc Mi ^ ^m) , i = A, . . . , D . (4.25) 



Hence, the wash-out processes are more efficient for larger neutrino masses, i.e. the final asymmetry 
should grow with the neutrino mass Mi. The finally generated asymmetry is not affected by the 
stronger wash-out processes, since for small rh\ the neutrinos decay late, where one can neglect the 
lepton number violating scatterings. 



This change of sign in the asymmetry is not observed in fig. |4.1| b. Due to the larger rh\ value 
the neutrinos are brought into equilibrium at much higher temperatures, where decays and inverse 
decays are suppressed by a time dilatation factor, i.e. the (s)neutrinos are produced in CP invariant 
scatterings off a (s)top. 



Chapter 5 

SO (10) Unification and Neutrino 
Mixing 



In order to study the implications of leptogenesis for low-energy neutrino physics and leptonic 
flavour mixing we will assume a similar pattern of masses and mixings for the leptons and the 



quarks in this chapter [15,17]. 



5.1 Neutrino Masses and Mixings 

If we choose a basis where the Major ana mass matrix M and the Dirac mass matrix mi for the 
charged leptons are diagonal with real and positive eigenvalues, 



(5.1) 
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the Dirac mass matrix of the neutrinos can be written in the form 



m T 



V 



\ 



m 2 
V m 3 / 



(5.2) 



where V and U are unitary matrices and the m; are real and positive. 

Since the Majorana masses M are assumed to be much larger than the Dirac masses m D , we 



have 6 Majorana neutrinos as mass eigenstates [11]. In the weak eigenstate basis the mass matrix 
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of the light neutrinos reads [JEJ] 

It can be diagonalized by a unitary matrix K, i.e. the light mass eigenstates 
have masses 



e,n,r , 



( 












m u 


V 









\ 


m VT J 



(5.3) 



(5.4) 



(5.5) 



whereas the heavy neutrino mass eigenstates 



l,R ' 



1,2,3, 



have masses 



M + O 



1 

17 



(5.6) 



(5.7) 



As we have seen in chapter ||, all the quantities relevant for baryogenesis, i.e. the decay widths, 
CP asymmetries and scattering cross sections, depend only on the product m\,m D , where the 
mixing matrix V drops out. On the other hand, the mixing matrix K in the leptonic charged 
current depends on the parameters of both unitary matrices U and V. Hence, leptonic mixing 
and CP violation at high and low energies are to leading order independent, i.e. the CP violation 
needed for baryogenesis does not allow to infer on CP violating interactions of light leptons. 

The mixing matrix U can be parametrized by three mixing angles and six phases. Five of these 
phases can be factored out with the Gell-Mann matrices Aj, 



U 



C7i e 



«A3<r „iA g r 



(5.8) 



In analogy to the Cabibbo-Kobayashi-Maskawa (CKM) matrix for quarks the remaining matrix U\ 
depends on three mixing angles and one phase. In unified theories based on SO(10) it is natural 
to assume a similar pattern of masses and mixings for leptons and quarks. This suggests the 
Wolfenstein parametrization [41] as an ansatz for U\, 

A\ z (p - irj) » 



/ 



Ui 



2 



-A 



A 



1 2 



y AX 3 (l -p-irj) -AX 2 



AX 2 
1 



(5.9) 



5.1. NEUTRINO MASSES AND MIXINGS 



59 



where A and \p + irj\ are of order one, while the mixing parameter A is assumed to be small. For 
the Dirac masses m$, SO(10) unification motivates a hierarchy like for the up-type quarks, 

mi = bX 4 m 3 ttt-2 = c\ 2 m 3 b,c = 0(l). (5.10) 



We have mentioned in section 4.1 that we also need a hierarchy in the Majorana masses Mj to get 
a non- vanishing lepton asymmetry. We choose a similar hierarchy as in eq. ( |5.10| ), 

Mi = BX 4 M 3 M 2 = CA 2 M 3 B,C = 0(1) . (5.11) 

Later on we will vary the parameters B and C to investigate different hierarchies for the right- 
handed neutrinos. 

Diagonalizing the neutrino mass matrix ( |5.3|) in powers of A yields the light neutrino masses 

b 2 

m,, = 



m v 



C + e^B\ X m ^ + 0(\*) , (5.12) 
A 2 rrv + (A 4 ) , (5-13) 



c 2 \C + e Aia B 



BC 

■2 



^ + 0(A 4 ) • (5.14) 



M 3 

We will not discuss the masses of the light scalar neutrinos here, since they depend on unknown 
soft breaking terms. 



In section 4.2 we have seen that the lepton asymmetry is largely determined by the mass 
parameter fhx, which is given by 

c 2 + A 2\ p + irj] 2 _ C ( c 2 + A 2 \p + ir l \ 2 ) 

B r ~~ c 2 \C + e Aia B\ " ' ' ' 

i.e. 777i is of the same order as the mass. According to eq. ( |3.39| ) the CP asymmetry in the 
decay of the lightest right-handed neutrino reads 



3 B A 2 



2 



4 777 3 

e l — ~ o To - ! ^To " — T -*- m 

'4 



8TT C 2 + A 2 \p + 77/| 2 v\ 



( p _ ir? )2 e ,2(«+V3/3)l + Q ( A 6) _ (5 16) 



In the next section we will always assume maximal phases, i.e. we will set 

_ 3 BA 2 lp + 77/l 2 4 7T7 2 6 
£ 1" - fr C 2 +A 2 \ p + lr ,\2 X ^2 +°{ X ) ■ ( 5 - 17 ) 

Hence, the lepton asymmetries that we are going to calculate may be viewed as upper bounds on 
the attainable asymmetries. 

Like in the non-supersymmetric scenario a large value of the Yukawa-coupling m^jvi will be 
preferred by this baryogenesis mechanism, since E\ oc m\jv\. This holds irrespective of our ansatz 
for neutrino mixings. 
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5.2 Numerical Results 



The neutrino masses ( 5,12 )-( ^.14j) can be used to constrain the free parameters of our ansatz. The 
strongest hint for a non-vanishing neutrino mass being the solar neutrino deficit^, we will fix the 



Vn mass to the value preferred by the Mikheyev-Smirnov-Wolfenstein (MSW) solution [43], 

m„ M ~ 3 • 10~ 3 eV . (5.18) 
Hence the parameter mj, which is of the same order as m v according to eq. ( p. 15 ), will be in the 



interval allowed by fig. |4.2 . 

The most obvious parameter choice is to take all 0(1) parameters equal to one and to fix A to 
a similar value as the A parameter of the quark mixing matrix, 



A = B = C = b = c=\p + ir]\~l 
A ~ 0.1 . 



The Ufj, mass in eqs. ( 5.13 ) and (|5.18|) then fixes the v e and v T masses, 



m v& a 8 • 10~ 6 eV , m Vr ~ 0.15 eV , 



(5.19) 
(5.20) 

(5.21) 



and m-i reads 



mi ~ 3 • 10~ 3 eV . (5.22) 
SO(10) unification suggests that the Dirac neutrino mass 777,3 is equal to the top-quark mass, 



7773 = m t ~ 174 GeV . 



This leads to a large Majorana mass scale for the right-handed neutrinos, 



M 3 ~ 2 • 10 14 GeV 



Mi ~ 2 ■ 10 10 GeV and M 2 ~ 2 • 10 1Z GeV , 



12 



(5.23) 



(5.24) 



and eq. (|5.17|) gives the CP asymmetry e\ ~ —6 -10 6 . Integration of the Boltzmann equations 
yields the ( B — L) asymmetry (cf. fig. |5.1| a) 



Y B - 



1 • 10" 



(5.25) 



which is of the correct order of magnitude. It is interesting to note that in the non-supersymmetric 
scenario one has Yb-l — 9 • 10~ 10 for the same choice of parameters. 



1 For a review and references, see Wm 
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Figure 5.1: Generated asymmetry if one assumes a similar pattern of masses and mixings for the 
leptons and the quarks. In both figures we have A = 0.1 and = rat (a>) and = nib (°)- 



Our assumption (5.23), m.3 ~ rrit led to a large Majorana mass scale M3 in eq. ( |5.24 ). To check 
the sensitivity of our result for the baryon asymmetry on this choice, we have envisaged a lower 
Dirac mass scale 

m 3 = m b ~ 4.5 GeV , (5.26) 



while keeping all other parameters in eqs. ( 5.19 ) and ( |5.20| ) fixed. The assumed mass ( |5,18 ) 
then yields a much lower value for the Majorana mass scale, 



M 3 ~ 1 • 10 11 GeV 



Mi ~ 1 • 10 7 GeV and M 2 ~ 1 • 10 9 GeV 



(5.27) 



while the light neutrino masses ( |5.21| ) remain unchanged. The CP asymmetry e\ — —4 • 10 9 
becomes very small. Consequently, the generated baryon asymmetry (cf. fig. |5,l|b) 



Y B -l ^1-10 



-12 



(5.28) 



is too small by two orders of magnitude. We can conclude that high values for both masses 
7713 and M3 are preferred. This suggests that {B — L) is already broken at the unification scale 
A GUT ~ 10 16 GeV, without any intermediate scale of symmetry breaking, which is natural in 
SO(10) unification. Alternatively, a Majorana mass scale of the order of 10 12 to 10 14 GeV can also 
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Figure 5.2: Solution of the Boltzmann equations suggested by the atmospheric neutrino problem. 

be generated radiatively if SO(10) is broken into SU(5) at some high scale between 10 16 GeV and 
the Planck-scale, and SU(5) is subsequently broken into the MSSM gauge group at the usual GUT 
scale ~ 10 16 GeV g§. 

Such a Majorana mass scale naturally leads to a baryogenesis scale M\ ~ 10 10 GeV. As discussed 



in section 4.2 and as one can see in fig. 5.1, the neutrinos can be brought into equilibrium at tem- 
peratures slightly above their mass, i.e. this scenario requires a reheating temperature ~ 10 10 GeV 
at the end of inflation. This is well compatible with the constraints on the reheating temperature 



from the gravitino problem [45| 



It is interesting to note that the v T mass in eq. ( |5.21| ) has got almost the value suggested by 
the atmospheric neutrino problem |}46|| 

m Vr tt7- 10" 2 eV , (5.29) 

if one assumes that the anomalous /i to e ratio produced by atmospheric neutrinos is due to 
oscillations from to v T , and when the mass is again given by the MSW value in eq. ( |5.18| ). 
If we then use eq. (HH) to fix the v T mass, the Dirac mass scale (HH) det ermines the Majorana 
mass scale 

M 3 ~ 4 • 10 14 GeV . (5.30) 
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Figure 5.3: Generated lepton asymmetry if one assumes a similar mass hierarchy for the right- 
handed neutrinos and the down-type quarks. 

The ratio of the v T and masses then yields the mixing parameter 

A = 0.15, (5.31) 
if the 0(1) parameters are given by eq. ( |5.19| ). The remaining neutrino masses read 

Mi ~ 2 ■ 10 11 GeV , M 2 ~ 9 • 10 12 GeV , m Ve ~ 2 • 1(T 5 eV . (5.32) 
Consequently, we again get a large CP asymmetry E\ ~ —3 • 10~ 5 , and a large [B — L) asymmetry 

(cf. fig. iiD 

Yb—l k 5 • 10~ 9 . (5.33) 

This shows that the parameters required to explain the anomalous results of neutrino experiments 
by neutrino oscillations also predict a baryon asymmetry of the correct order of magnitude, although 
the large mixing angle which seems to be required to solve the atmospheric neutrino problem is 
difficult to accommodate within our small-mixing ansatz. 

Up to now we have always assumed a mass hierarchy for the heavy Majorana neutrinos like 
for the up-type quarks. Alternatively, one can assume a weaker hierarchy, like for the down-type 
quarks by choosing 

B = 10 and C = 3 . (5.34) 



64 



CHAPTER 5. SO (10) UNIFICATION AND NEUTRINO MIXING 



Keeping the other parameters in eqs. (5.19) and ( 5.20| ) unchanged fixes the v e and v T masses, 



m v ~ 5 ■ 10 _b eV 



m Ur ~ 0.7 eV , 



and the mass parameter mi, 



(5.35) 



(5.36) 



mi ~ 1 • 10 3 eV . 

Choosing the Dirac mass scale fl5j|D we get a large Major ana mass scale 

M 3 ~ 4 • 10 13 GeV Mi ~ 4 • 10 10 GeV and M 2 ~ 10 12 GeV . (5.37) 

From eq. ( |5,17| ) one obtains the CP asymmetry e\ — —6 • 10 -5 . The corresponding solutions of the 



Boltzmann equations are shown in fig. 5.3. The final (B — L) asymmetry. 



Yb—l ^ 2 ■ 10" 



(5.38) 



is much larger than requested, but this value can always be lowered by adjusting the unknown 
phases. Hence, the possibility to generate a lepton asymmetry does not depend on the special form 
of the mass hierarchy assumed for the right-handed neutrinos. 

In the non-supersymmetric scenario one finds for the same parameter choice 



Yb—l ^2-10- 



(5.39) 



Hence, when comparing the supersymmetric and the non-supersymmetric scenario, one sees that 
the larger CP asymmetry in the former and the additional contributions from the sneutrino decays 
are compensated by the wash-out processes which are stronger than in the latter. 



Conclusions 



We have analysed in detail the generation of a cosmological baryon asymmetry by out-of-equilibrium 
decays of heavy right-handed Majorana neutrinos and their scalar partners in a super symmetric 
extension of the Standard Model. By developing a resummed perturbative expansion in flavour 
non-diagonal self-energies, we could show how self-energy contributions to CP asymmetries in 
heavy neutrino decays are consistently taken into account. 

We have discussed all the decays and scattering processes relevant for leptogenesis, and by 
solving the Boltzmann equations we have shown that, in order to be consistent, one has to pay 
attention to two phenomena which can hamper the generation of a lepton asymmetry. 

First, one has to take into consideration lepton number violating scatterings mediated by a 
heavy neutrino or its scalar partner. These processes, which are usually neglected, can wash out 
a large part of the asymmetry if the Yukawa couplings of the right-handed neutrinos become too 
large. 

On the other hand, the neutrinos have to be brought into thermal equilibrium at high temper- 
atures. We could show that for this purpose it is not necessary to assume additional interactions of 
the right-handed neutrinos in our theory, since the Yukawa interactions can be sufficiently strong 
to produce a thermal population of heavy neutrinos at high temperatures, while still being weak 
enough to prevent the final asymmetry from being washed out. 

The observed baryon asymmetry can be obtained without any fine tuning of parameters if one 
assumes a similar pattern of mixings and Dirac masses for the neutrinos and the up-type quarks. 
Then the generated asymmetry is related to the mass, and fixing this mass to the value preferred 
by the MSW-solution to the solar neutrino problem leads to a baryon asymmetry of the requested 
order, provided (B — L) is broken at the unification scale, as suggested by supersymmetric SO (10) 
unification. The baryon asymmetry is generated at a scale of approximately 10 10 GeV, which looks 
promising with respect to the gravitino problem. 
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CONCLUSIONS 



In supersymmetric theories there are further possible sources of a (B — L) asymmetry, e.g. it 
may be possible to combine inflation with leptogenesis by using a right-handed scalar neutrino as 
the inflaton (cf. refs. [^Tj). In this connection, possible constraints on the neutrino masses and 
on the reheating temperature from lepton number violating processes at low temperatures require 
further studies. 

Furthermore, it should be studied to which extent a "primordial" baryon asymmetry, i.e. an 
asymmetry generated during or shortly after reheating, is affected if right-handed neutrinos come 
into equilibrium after reheating. This may yield interesting constraints on Yukawa interactions of 
first generation leptons, since an asymmetry in right-handed electrons might be protected from 
being washed-out if right-handed electrons decouple from the thermal plasma 

Finally, as discussed in the first chapter, one has to separate "on-shell" and "off-shell" contri- 
butions to lepton number violating scatterings, in order to be able to describe the generation of 
a baryon asymmetry in terms of Boltzmann equations. To avoid this separation and treat these 
contributions simultaneously, one has to go beyond the semi-classical approximation realized in the 
Boltzmann equations by constructing a complete nonequilibrium quantum kinetic theory. Although 
such a theory has not been realized up to now, quantum corrections to the Boltzmann equations 



have recently been investigated in refs. 4£]. 



Appendix A 

One-Loop Integrals 



We summarize some standard formulae for dimensionally regularized one-loop integrals in 
Minkowski space. We follow the notation of refs. |5(J, although we use a metric g^ u = (+, — , — , — ) 
(cf. e.g. ref. §1). 



A.l One-Point Function 



In n = 4 — 2e dimensions the scalar one-point function is defined by 



u A ~ n r 1 
A(mi) = L / d n k T -— . (A.l) 

itt z J k l — m{ + is 



In the limit e — > 0, A{m\) is given by 

A( mi ) =ml(A- In + , (A.2) 

where the UV-divergence is contained in 

A = - - C + ln(4vr) , (A.3) 

and C = 0.577216 is Euler's constant. Note that the massless tadpole ^4(0) vanishes in dimensional 
regularization, and that A{m\) has no absorptive contribution 

Im [A(mi)\ = . (A.4) 



A.2 Two-Point Functions 



Three different two-point integrals can occur 

4-n 

2 \ P 



B ( Pl , mi ,m 2 ) = j ^ {k2 _ ml+ie)[{k+pi)2 _ ml + i£] (A.5) 
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B fJ ,(pl,m 1 ,m 2 ) 

, mi, 7712 J 



4— n 



77T^ 
4— n 



77T^ 



(A; 2 — m 2 + ie)[(k + pi) 2 — m 2 + ie] 
kfj,ki/ 

(k 2 — m 2 + 7e)[(/c + pi) 2 — ml + 



(A.6) 
(A.7) 



Lorentz covariance of the integrals allows to decompose the tensor integrals into tensors constructed 
from the external momentum p\ , and the metric tensor g^ v 

B ll (pj,m 1 ,m 2 ) = pi^B 1 (pl,m 1 ,m 2 ) , (A. 8) 

Bfj,u(pl,mi,m 2 ) = Pi,nPi,vB 2 i(p(, mi, m 2 ) + g^ u B 22 (p{, mi, m 2 ) . (A.9) 

Using the Feynman parametrization, one can derive an integral representation for Bq, 



Bo(p\,rrii, m 2 ) = A — J dx In 
b 



9 9 / 9 9 9\ o 

x Pi — x(pi + mf — m^j + — is 



+ 0(n - 4) , (A.10) 



which yields the following useful identities in the limit ti — > 4 
B (p?,0,0) = A - In (&j + 2 + iir6(p 2 ) 



B (0,0,m) = B (0,m,0) = A-ln(^) + 1 = -^(m 2 ) . 

V A* / m z 



(a.ii; 



(A.12) 



Contracting eqs. ([A.6|) -( [A9| ) with pi )(U and yields a set of coupled linear equations, which 
determine the scalar coefficients B\, B 2 \ and B 22 , 



£i(Pi,mi,m 2 ; 



S 2 i(^,mi,m 2 ; 



S 22 (p 2 ,mi,m 2 ; 



— ^ \A{mx) - A{m 2 ) + (m 2 2 - m\ - p 2 )B (p 2 , mi, m 2 ] 
2pf 



1 

~2~e 

1 
M 



+ UV-finite parts, 



■2 [A(m 2 ) - mlB (pl,mi,m 2 ) 
-2{p\ + m\- ml)B 1 (p 2 ,m 1 ,m 2 ) - -{m\ + m 2 . - ~pfj 

h UV-finite parts, 

3e 



- [A(m 2 ) + 2mjBo(p 2 ,mi,m 2 ) 
+(p\ +m\ - m 2 )Bi(p\,mi,m 2 ) + m 2 + m 2 . - -p\ 



12e 



(p 2 — 3m 2 — 37712) + UV-finite parts. 



(A.13) 
(A.14) 

(A.15) 
(A.16) 

(A.17) 
(A.18) 
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In the case of equal or vanishing masses one has 



Bi(pj,m,m) = --B (pl,m,m) , 



B 2 i(p?,0,0) = - 



B (plAO) + - 



-Pi 
12 



5 2 2(P?,0,0) 



A. 3 Three-Point Functions 



In general one has four different three-point functions 



6 



Bo(pl,0,0) + - 



C (pl,pl,m 1 ,m 2 ,m 3 ) 



4— n 



ITT 
4— n 



s 



C fl {p 1 ,p 2 ,m l ,m 2 ,m 3 ) = 



Cn„{pl,pl,m 1 ,m2,m 3 ) 
Cn„ p (pl,P2,m 1 ,m2,m 3 ) 
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4—n 



d n k 
d n k 
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Dm\ Dm2 (Pi ) 


(Pi 


P2) 


kp 
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P2) 


k/j, k v 






D mi Dm 2 {Pi)D m:i 


(Pi 


P2) 


k^kykp 






D mi D m2 (pi)D m3 


(Pi 


P2) 



where we have introduced the following abbreviations for inverse propagators 

D mi = k 2 - m\ + ie , 
D m2 (Pi) = (k + pi) 2 - ml + is , 
D m Jpi,P2) = (k+p 1 +p 2 ) 2 -ml + ie . 



(A.19) 
(A.20) 

(A.21) 



(A.22) 
(A.23) 
(A.24) 
(A.25) 

(A.26) 
(A.27) 
(A.28) 



Like for the two point-functions, Lorentz covariance of the integrals suggests the following tensor 
decomposition of the tensor three-point functions 



Cn(pl,pl,m 1 ,m 2 ,m 3 ) = pi^Cn + P2,pC vl , 

C fll/ {pl,pl,m 1 ,m 2 ,m 3 ) = pi,pPi,„C 21 + p 2: pp 2>u C 22 

+ (PlP2)( MJ /)C23 + 9^C 2 4 , 
C fll/p {pl,pl,m 1 ,m 2 ,m 3 ) = pi,pPi,uPi, P C 31 + p 2: pp 2:l/ p 2:P C 32 

+ {P2P\Pl)(pup)C 33 + (PlP2P2)(pup)C 3 4 
+ (Pig)(pvp)C35 + {P29)(pup)C 3 Q , 



(A.29) 



(A.30) 



(A.31) 
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where we have used the following abbreviations for index symmetrizations 

(PlP2)( M u) = Pl,pP2,u +PX,vP2,n ! 
(PlP2P2)( M up) = Pl,i*P2,vP2,p + P2,»Pl,vP2,p + P2,iiP2,uPl,p , 
(pig)(pup) = Pi,p9v P + Pi,v9p,p + Pi, P 9fiv ■ 



(A.32) 
(A.33) 
(A.34) 



The form factors can be related to the scalar functions A, Bq and Cq by contracting the 
definitions ( A.23| )-( |A.25 ) and ( A.29| )-( [A.31j ) with external momenta p\, P2 and the metric g^ v . 

In our calculation we only need three-point functions with two vanishing masses (m 2 = ITI3 = 0), 
and two light-like momenta p\ = and {p\ + P2) 2 = 0. Then C12 an< ^ C11 re& d 



Cu(pi,P2,m,0,0) 



1 



2pi ■ p 2 

C u {pi,P2,m,0,0) = 2Ci 2 (pi,P2,m,0,0) , 



[B (0,m,0) - B (pl0,0) -m 2 Co{ Pl ,p2,m,0,0)] , (A.35) 

(A.36) 



and the imaginary parts of Cq and C12 are given by 
Im[C (pi,P2,m,0,0)] : 

Im[Ci 2 (pi,P2,m,0,0)] : 





-In 


fi + ii 


Pi 




\ m 2 




1 - 


m 2 / 


p\ 







pi 

m 2 



(A.37) 
(A.38) 



Appendix B 

Spinor Notation and Conventions 



We will use the conventions of ref. [30] with flat space-time metric = (+,—,—,—). Greek indices 
a, (3, a and j3 run from one to two and denote two-component Weyl spinors, while all other Greek 
letters denote Lorentz-indices. 

B.l Weyl Spinors 

Two-component spinors ip and tp transform under the (^,0) and (0,^) representations of the 
Lorentz group SO(l,3). Matrix representations are given by the universal covering group SL(2,C) 
of SO(l,3), i.e. under a Lorentz transformation M £ SL(2,C) the Weyl spinors transform like 

r// a = M-^ V , = (M*)- 1 ■ " ^ ■ (B.2) 

The Pauli matrices o ^ form a basis of SL(2,C), 

Spinor indices can be raised and lowered using the antisymmetric e-tensors 

- = i!o)' ^=(- iii- (R5) 
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a/3 

i.e. one has 



10/ \ -1 



^ = £ a % , V'a = e aP ^ , (B.7) 

^ = e^iP- , Y>- = . (B.8) 

Correspondingly, spinor indices of Pauli matrices can be raised and lowered, 

a^ a = e^e^a^ . , (B.9) 

i.e. the matrices cf M are given by 

a° = a° , a 1 ' 2 ' 3 = -a 1 ' 2 ' 3 . (B.10) 



The Pauli matrices fulfil 



(a^a u + a u a»)f = 2g^Sj , (B.ll) 

{a»a v + = 2g> u, S Sl ., (B.12) 

Tra'V = 2(^", (B.13) 

= 2(5/5/. (B.14) 
Products of two-component spinors are defined such that 

^x = 4> a x a = -i>ax a = x a i> a = X4> , (B.15) 

tpX = = -^"x^ = X^ a = XV' , (B.16) 

(x# = (x Q VO f ^X^^X- (B.17) 
Products of Weyl-spinors involving Pauli matrices read 

XCJ ^ = -Y^x, (B.18) 

(x^?) f = V"^X, (B.19) 

X^cf^ = i/h^Wx, (B.20) 

(X^a» f = ^VL (B.21) 

(^A)x- = ^(A^x)W^)a • (B- 22 ) 



B.2. FOUR-COMPONENT SPINORS 
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Furthermore, when computing superfield products (cf. app. B.3), one can take advantage of the 
relations 

(B.23) 
(B.24) 
(B.25) 
(B.26) 
(B.27) 

B.2 Four- Component Spinors 

We define the Dirac 7-matrices as 

-\ ° a ] , 7 5 =WtV= I _1 ° I (B.2-S) 
/ V01 



aaa/3 _ 

V 1/ — 




- F a/3 2 
c u 






2 


6 a 6f3 = 


1 

2 


e a/ 3 9 2 , 


Te h = 


1 






2 










e-e- = 






a p 




2 «/3 


a^ee^e = 


1 

2 





obeying 

{7^,7-} = 2«T, (B.29) 

{7 5 ,7 At } = 0, (B.30) 

(7 5 ) 2 = 1 • (B.31) 

As usual, 7 intertwines the 7^ representation of the Dirac algebra with the equivalent hermitian 
conjugated representation 7^ 

7°7 M 7° = 7 Mt , (B.32) 

7°7 5 7° = - 7 5t = -7 5 . (B.33) 

This representation of the 7-matrices can be used to relate Weyl spinors to the more familiar 
four-component spinors. A Dirac spinor ^ D consists of two Weyl spinors 

i.e. its hermitian conjugate reads 

*d = ^7°= (A", X-) • (B.35) 
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The chiral projectors read 




1 



1 




(B.36) 



i.e. for the Dirac spinor (B.34) one has 



The charge conjugation matrix C is defined by 



2„,0 



C = —i^f 7' 



and it fulfils the following useful identities 



£ a f3 

e¥ 



c T = c ] = c- 1 = -C , C 2 



-1 . 



C intertwines the 7^ and — 7J representations of the Dirac algebra 



C^C- 1 = 7 5 T 



(B.37) 
(B.38) 



(B.39) 



(B.40) 



(B.41) 
(B.42) 



The charge conjugated Dirac spinor (B.34) then reads 



A a 



i.e. one has 



On the other hand, a Majorana spinor ^ M contains only one Weyl spinor 



4% 



i.e. eq. ( |B.43| ) immediately implies 



(B.43) 



(B.44) 



(B.45) 



(B.46) 



B.3. SUPERFIELD PRODUCTS 
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In Lagrange densities we can switch from Weyl spinors to four-component spinors and vice versa 
by means of the following relations 



~^P L ^ 2 = Mx 2 , 



(B.47) 
(B.48) 
(B.49) 
(B.50) 
(B.51) 
(B.52) 



where Vl/j (i = 1, 2) is a generic Dirac spinor consisting of two Weyl spinors x- an d \ (cf- e Q- ( B.34D ). 



B.3 Superfield Products 

Since products of chiral superfields are again chiral, these products can be conveniently computed 
in the y-basis ( p. 31 ) . For the products of two or three chiral superfields <!> i one finds 



$i(y)$ i (y) = Mv)Mv) + v^fl Mv)My) + MWM 



+e 2 



My)Fj(y) + My) F i(y) - A(y)i>j(y) 



$i(.y)$j(y)$k(y) = A i (y)A j (y)A k {y) 



(B.53) 



(B.54) 



+V26 



Mv)My)Mv) + Mv)Mv)Mv)+1>k(v)Mv)Mv) 



+9 2 F i (y)A j (y)A k (y) + F^A^A^y) + F k {y)Ai(y)Aj(y) 



-A(y)' i Pj(y)Ak{y) - i>j{y)^k{y)Ai{y) - i>k(y)i>i(y) A j(y) 



On the other hand, the product of a chiral superfield &j and an antichiral superfield $| is neither 
chiral nor antichiral. In terms of the variables x^, 6 and 9 it reads 



$i{x)$j(x) = A*{x)A j (x) + V^ieiPjix^AUx) + y/2(ptl)i{x))Aj{x) 



(B.55) 



e 2 A*{ x )F 3 (x) + e*A 3 {x)F;{x) 



ia a ^ (A*(x)d ll A J (x) - A J (x)d^A*(x)) - 2^(sty it 
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+e 2 e' 



772 



+e 2 0- 



-jLaJ* (A*(x)d^?(x) - ^(x)d,A*(x)) - V2F 3 {x)^{x) 



-jfaa" Ui a (x)d ti A j (x) - A J (x)d^ l a (x)j + V2F* (x)^ Q (x) 
F*(x)F J (x) - ±A*(x)nMx) - - A A 3 {x)uA*(x) 



+\{d tl A* i {x)){d»A j {x)) - l-d^xp^x) + % -Ux)a^d^ 3 {x) 



Superspace integration discussed in Chapter ^ will project out F-terms, i.e. terms proportional to 2 



in chiral superfields and the D-term (proportional to 9 2 6 ) in eq. ( B.55 ). These are the terms which 
can be used to construct SUSY-invariant actions since they transform into a spacetime derivative 
under SUSY transformations. 

When computing 2^2 scatterings in superspace one has to evaluate products of four chiral 
or anti-chiral superfields, with or without covariant derivatives acting on them. In our calculations 



we need the following two products, which can be computed by successively using eqs. ( B.53 ) and 

(HI) 



f d 2 e d 2 e QrixuO, )*,(*!, 0, 0) X -~ e, e ) $,(x 2 , o,o) = 

= A r (x l )F s {xi) + A s (x l )F r (xi) - ip r (x 1 )Tp s (xi) A i (x 2 )A j (x 2 ) 
+ Ai(x 2 )Fj(x 2 ) + Aj(x 2 )Fi(x 2 ) - ipi(x 2 )ipj(x 2 ) A r (x 1 )A s (xi) 
+ i) s {xi)A r (x\) +ip r (xi)A a (xi) ifji(x 2 )Aj(x 2 ) +ipj(x 2 )Ai(x 2 ) 



(B.56) 



J d 2 e d 2 e $ r (xi,e, ) 0, ) $i(x 2 , 0, e ) $j(x 2 ,e, ) 



-A\(x 2 )A\{x 2 )UiA r {xi)A s {xx) 



+ 1 
+ 



Wi(x 2 )A]{x 2 ) +Tj(x2)Al(x 2 ) 
A r {xi)F s (x l ) + A s (x 1 )F r {xi) - ip r (x 1 )ip s (x 1 ) 
A){x 2 )F}{x 2 ) + A\{x 2 )f]{x 2 ) -Tj(x 2 )Wi(x2) 



(B.57) 



x 



■^ s (x 1 )A r {xi) + ip r (x 1 )A s (x 1 ) 



Here we have partially integrated derivatives and dropped total derivatives which do not contribute 
to the action. 



Appendix C 

Feynman Rules 



In this appendix we present the component field Feynman rules that we have used in the calculations 
of chapters [l] and ^. Feynman rules for Majorana fermions which yield the correct relative minus 
signs between different diagrams contributing to a process without explicit recourse to Wick's 



theorem are proposed in ref. [52|. The basic idea is to introduce a continuous fermion flow, i.e. 
an arbitrary orientation of each fermion line. Then one can form chains of Dirac matrices by 
proceeding in a direction opposite to the chosen fermion flow. Relative signs of interfering diagrams 
are determined like for Dirac fermions, i.e. any permutation of two external fermion lines gives a 
minus sign. One only needs the Dirac propagator for all fermions. However, one has to introduce 
two analytical expressions for each vertex involving fermions, corresponding to the two different 
choices of the fermion flow. 

As an example, consider the coupling of a right-handed neutrino to a SM Higgs and lepton 
doublet. This coupling can be written in two equivalent ways 

N\ T v P^(leH 2 ^ =- (H 2 el^P h X v N , (C.l) 

corresponding to the two possible choices for the fermion flow. This gives rise to the following 
equivalent vertices, where the thin arrow denotes the chosen fermion flow 

I I 

y y 
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In general, one has spinor structures of the kind 



(C.2) 



where the four-component spinors tpi and ip2 are either Dirac or Majorana fermions, and T is a 
product of Dirac 7-matrices, 

r = i,y\ 7 5 )7 y>^- (c.3) 

Reverting the fermion flow corresponds to replacing particles by antiparticles, i.e. ( |C.2| ) is rewritten 
in the equivalent form 

where 

r for r = i,7 5 , 7 V 

-r for r = 7^,0-^ 

Hence, when stating Feynman rules in the following we can restrict ourselves to one fermion flow. 
Changing the fermion flow just amounts to replacing V by ±r, according to eq. ( |C.5| ). 

Decomposing the superfield products in the superpotential ( [O] ) into component fields, we get 
the Yukawa interactions of a right-handed Majorana neutrino 



T' = CT T C~ l 



(C.4) 



(C.5) 




-{iK) l3 e ah (Pj Po 




1. a 



Correspondingly, the interactions of a scalar neutrino are given by 



1 

y 



(iX^ijS^iP^p 



N? 
1 

A 
1 
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Here we have not specified an explicit fermion flow, since these diagrams have a natural orientation 
of the fermion lines. 

The mass term in the auxiliary neutrino field F N c (cf. eq. ( 3.22| )) yields trilinear scalar couplings 



i 

* 
i 

A 

1 A 



-{i\ v M)ij e 



ab 



k 
i 

A 
i 

A 

y a 



-(iMXDjie' 



t \ . . ,06 



Furthermore, we take into account the following Yukawa couplings of the (s)top 



7 b 

H,f3 



I 

* 
I 




UJ 
I 

A 
I 






„b 

H,/3 



i 

y 




-(a u )ye a6 (P L )a/3 



jra 
I 

A 
I 




-(iXl^e^iP^c 



Finally, the scalar potential 



V = Y J F l F i, i = H u H 2 ,Q i ,...,N i 



c 

i > 



(C.6) 



yields quartic scalar couplings involving one or two scalar neutrinos 



< y 



- a 



\ / 1 

< y 



-i(At)ji(A u ) sr 5 a,) 
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jja 



\ ✓ 

' S v 



*(a+a,)^ 



3 

jja 



\ ✓ 
\/ 
✓ \ 

J * 



-i{\t\ v )ij 5 ab 



S a / 3 

< -f 
\ / 

✓ \ 

' S v 

\ 

E c ' s Ho 

-b 



i(XtXi)jid ab 



~a 

L \ 



' I, 



\/ 

J i 



—i{\i/) r j (Xl)is <5 ab 



Internal Dirac or Majorana fermion lines are all represented by the usual Dirac propagator 

P f i 



a 







J-m + ieJ p a 

Correspondingly, we assign spinors to external fermion lines with orientation 

• > ; < ; • «(P, s) , 

• ; * • ; * • ; v (p, s ) , 

» < ; ' ; 1 u (p> s ) > 

; > ; < ; *>(p> s ) , 

where the momentum p always flows from left to right. The spinors v and u are related by charge 
conjugation 



v(p, s) = C u T (p, s) 44> u(p, s) = v T (p, s) C 
u(p, s) = Cv T (p, s) 44> v(p, s) = u T (p, s) C 



(C.7) 
(C.8) 



Appendix D 



Kinetic Theory 



The microscopic evolution of particle densities and asymmetries is governed by a network of Boltz- 
mann equations. In the following we will compile some basic formulae to introduce our notation 



D.l Thermodynamics in the Expanding Universe 

The early universe can be assumed to be spatially homogeneous and isotropic. Hence, it is described 
by a Robertson- Walker metric 

ds 2 = dt 2 - R(t) 2 / dr j _ + r 2 d9 2 + r 2 sin 2 9d6 2 \ , (D.l) 
[ 1 — kr z J 

where (t,r,6,(f)) are comoving coordinates and k = ±1,0 describes the spatial curvature of space- 
time. The scale factor R(t), which describes the expansion of the universe, is given by the Friedmann 
equation 

R +k = ^pR" , (D.2) 
where p is the energy density of the universe, and 

G=^r (D.3) 
denotes Newton's constant in units where H = c = 1, and m Pt = 1.2211-10 19 GeV is the Planck mass. 



Neglecting the curvature term k in the Friedmann equation (D.2), which is a good approximation 



in the early universe, we get an equation for the Hubble parameter H 



H = — - -LJ*2. . (D.4) 
R m Pl V 3 v ' 
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In a radiation dominated universe the energy density reads 



P = g*-T 4 , (D.5) 



where g* is the number of effectively massless degrees of freedom 

g*= Yl 9i + s ^ 9i ' ( D ' 6 - ) 

i— Bosons i — Fcrmions 

and gi is the number of internal degrees of freeedom of the corresponding particle. At temperatures 
far above the electroweak scale one has g* = 106.75 in the standard model, and g* = 228.75 in the 
MSSM. 

Hence, the Hubble parameter in a radiation dominated universe reads 



D.2 Boltzmann Equations 

It is usually a good approximation to assume Maxwell-Boltzmann statistics, so that the equilibrium 
number density of a particle i is given by 

n^T) = j£^J d 3 K r with fP (E h T) = e~^ T . (D.8) 

For a massive particle one finds 

„ r( T) = ^ (! (2) . (D.9) 
whereas for a massless particle one gets 

nf(T) = ^ . (D.10) 

Particle densities can be changed by interactions and by the expansion of the universe. Since 
we are only interested in the effect of interactions, it is useful to scale out the expansion. This is 
done by taking the number of particles per comoving volume element, i.e. the ratio of the particle 
density rii to the entropy density s, 

Y t = ^, (D.ll) 

s 

as independent variable instead of the number density. In a radiation dominated universe the 
entropy density reads 

2tt 2 

s = g*—T\ (D.12) 



D.2. BOLTZMANN EQUATIONS 
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In our case elastic scatterings, which can only change the phase space distributions but not the 
particle densities, occur at a much higher rate than inelastic processes. Therefore, we can assume 
kinetic equilibrium, so that the phase space densities are given by 



fi(Ei,T) 



eq 



-Ei/T 



(D.13) 



In this framework the Boltzmann equation describing the evolution of a particle number in an 



isentropically expanding universe reads [36,13] 



d^ 
dz 



sH (m^) 



E 

a,i,j,.. 



Y^Yg . . . 



YY 



y-eqx^eq 
1 i 1 -j ••• 



ip + a + ...) 



(D.14) 



where z = m^/T and H (m^) is the Hubble parameter at T = m,^. The 7 eq are space time densities 
of scatterings for the different processes. For a decay one finds Jl3| 

(D.15) 



7D:=7 eq (V<-i + j + ...)=< qKl(Z) 



* K 2 (z) ' 

where Ki and K2 are modified Bessel functions and T is the usual decay width in the rest system 
of the decaying particle. Neglecting a possible CP violation, one finds the same reaction density 
for the inverse decay. 

The reaction density for a two body scattering reads 

00 _ 

^q (V , + a ^ + J - + ...) = J_ J ds t( s) ^ Kl (^pj , (D.16) 

where s is the squared centre of mass energy and the reduced cross section a(s) for the process 
ip + a — *i + j + ... is related to the usual total cross section a(s) by 

2X(s,ml,m 2 a 



a(s) , 



(D.17) 



where A is the usual kinematical function 



A(s, rri^,ml ) = [s - (m^ + m a ) 2 ] [s - (m^ - ma) 2 ] 



(D.18) 



Appendix E 



Reduced Cross Sections 



In this section we will collect the reduced cross sections for all the 2 «-> 2 and 2^3 processes 
that we had discussed in chapter |3[ The corresponding reaction densities, which can be calculated 
analytically in some interesting limiting cases, will be discussed in the next appendix. 

E.l Lepton Number Violating Processes Mediated by Right- 
Handed Neutrinos 



We have mentioned in the main text that we have to subtract the contributions coming from 
on-shell (s)neutrinos, i.e. we have to replace the usual propagators by off-shell propagators 

1 X — (In ,1 x — do 



Dj(x) ' (x-ajf + ajCj 



and 



Dj{x) ' {x - aj) 2 + ajcj 



(E.l) 



To begin with, let us specify the reduced cross sections for the reactions depicted in fig. 3.5. For 
the processes I + h *-* I + h and I + H2 «-> I + H\ one has 



X 

— V 



X 



+ 



X 



+ E Re 



n,3 



II J. 



+ (x + aj) 



^Ja n aj 



1 



+ 



aj Dj{x) 2D](x) 
x 2 

+ 



1 + 



x + aj . , 

dm 1 1,1 



x + a 



a j 



Dj(x) D n (x) Dj(x)D n (x) 
x + aj \ , , / 2 



(E.2) 



an aj i^ n y^ j j \ Li j j X^J -^j \ 

where n and j are the flavour indices of the neutrinos in the intermediate state. The interference 
terms with n 7^ j are always very small and can safely be neglected. 



D n {x) 



In 



- J + {x + a r< 



1 



Dj{x) 



In 



x + a r 
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The reduced cross section for the process I + h I + h\ reads 



+ 



x + o, £)j(x) 2£>J(x; 



+ 1 



D,(x) 



In 



x + a 



a,; 



sj Oj n CLj 



X X 
+ 



Dj(x) D n {x) Dj(x)D n (x) 



(E.3) 



2 1 

+ 



a n - a,j D n (x) 



In 



x + a. 



2 1 

+ 



n \a,j - a n Dj(x) 



In 



x + aj- 



ar. 



The same result is valid for the CP conjugated process. 

~ ~\ t 
For the process I + h <-+ I + B.\ one finds 



+E Re [( A ^)l- 



X 



X X 

+ — In 



aj (x + a 3 ) D ^ x) D .( 



x + a 



a; 



n,3 



2x 2 



D j (x)D n (x) 



+ x 



2 1 

+ 



In 



(In, CL i 



3 D n {x) 



x + a. 



+ x 



2 1 
+ 



In 



x + a r , 



a,, 



(E.4) 



For the scattering I + Hi — ► / + f/ c + g. and the corresponding CP transformed process we have 



+ 



J>e (At A,)'. 



XX X 

+ 



X + GU 
1+ J 



^(x) (x) V ^(*). 



In 



x + a. 



Q n CLj 



X 



X X 

+ — — + 



X 



Dj(x) D n (x) Dj{x)D n (x) 



(E.5) 



+ (X + Oj) 



In 



X + Oj 



+ (x + a T . 



In 



x + a r , 



a,, 



Finally, we have two processes which do not violate lepton number but merely transform leptons 
into scalar leptons and vice versa. We have the 2^2 scattering I + H2 I + h, 



3,n 



it J 



and the 2^3 process l + h^l + tf + U C , 

3a u 



£»,(x) J D n (x) ' 



2 



Dj(x)D n (x) 



(E.6) 



(E.7) 
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Let us now come to the 2^3 processes shown in fig. p.6| . 
For the transition q + U c — » / + I + H2 the reduced cross section reads 



a ; 



jj x 



■ I ^ ■( ill 



dj + Cj 



djCj 



arctan 



x — di 



V 



+ arctan ( . / -±- 



djCj 



In 



[X — dj) + djCj 
d) + djCj 



- I dxi — In 

2 I Dj( Xl ) 



(x — X\ — dj) 2 + djCj 
d 2 + djCj 



+ 2^Re 



OjfiQij 



- [ dx\ -=r- In 

2 J D n (x x ) 



(x — X\ — dj) 2 + djCj 



dj + djCj 



(E.8) 



+ 2 J djCj 



X — d r 



(dj - d n y 



/ X — dj 

arctan 1 



+ arctan . / 



X ~ a 3 ln [ {x-dj) 2 +djCj 



dj + djCj 



+ 27 



X — <2„ 



(d n - dj) 



• X Ob'i b , 

arctan | — — | + arctan 



X Q-n , / (x d n ) + d n C n 

H in 



d n -\- d n C r 



The remaining integral cannot be solved exactly. However, it can be neglected for x > dj,d n and 
for x < dj,d n it can be approximated by 



dxi 



D n (xi) 



■ln 



(x — X\ — dj) + djCj 
d 2 + djCj 



(E.9) 



ln 



ln 



dn 



+ Sp 



Sp 



(X n x 



where Sp(x) is the Spence function or dilogarithm 



Sp(x) = Li 2 (x) = - dy 



ln(l - y) 



(E.10) 



For the scatterings q + U c — > I + 1 + h, V + q ^ I + + h and V + U c — > I + cfi + h the reduced 
cross sections are equal, 



^(x)=^\x) = l^-{^(xiX 



8vr 2 : 



--x + -(x - 2oj)ln i -j 



(x — Oj) 2 + OjCj 



1 fa 



a V c 



X — dj 

arctan 1 



y/djCj 



+ arctan 



(E.ll) 
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+ 2 E 



n-,3 



a£a„ 



»j 



-2x + Oj 



x — a,- 



In 



(x — a,-) + a,c 



m 1 



+ 2w a,-c 



xa n - 2a„a-,- + a 2 



arctan 



x — a,- 



7 



+ arctan 



CljCj 



+ 2\fa r ^ n 



xa,j — 2a n a,j + a n 
(a„ - cij) 2 



x — a n \ I cl„ 

arctan ( — I + arctan J — 
Va n c n J \\ c n 



~t ~ ~ — t 

For the process I + q —* / + C/ c + H2 an d similar reactions one gets 



4 10) (.) 



3^ 
16vr 2 



"V V hi x 



In 



[x — a,) + a,c 



+ 



X — CLj 



V a i c i 
+ 2^Re 



arctan 



x — a, 



a£a„ 



v \/ a j c i 



+ arctan . / 



+ 2 / dxi 



£>j(si) 



2 a: + Oj ^ / x + Oj \ | 2 a: + a w ln ^ a; + a n 



Sp(--]-Spl-^ 



(2 jf d n 



(E.12) 



+ 



x — au 



In 



(x — a,,-) 2 + ajCj 



+ 2\/ ajCj- 



arctan 



X — CLj 



V a i c i 



+ arctan 



H in 1 



a; 

+ y 





O 2 + ClnCn 



+ 2^ 



x — a,- 



drt.Cn 



(a„ - a,-) 2 



arctan ( — | + arctan 



£>j(xi) V V a « 



s P (-^Y)+,J— 

«n// D n (xi) 



\/(l n C ri 



Sp j -— ) - Sp 1 



The remaining integral can again not be solved exactly. However, it can be approximated by 



dxi 



Si>!-- -Sp 



(E.13) 



x 



arctan 



x — a,- 



7 



+ arctan 



2a n \ aj + ajCj 



for x < a n and for x > a n it can be neglected. 



Finally, we have to compute the t- and u-channel processes in fig. 3.7, which give simple con- 
tributions. 
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For the processes I + I <-> h + h and I + I <-> H\ + h\ we get 



*T\ X ) = ^ 3) (-) = ^|E(^) 

+ E Re (xUv)]. 



x 



+ 



ttj In 



+ 



1 



x + aj x + 2a,- 



+ 



x + a 



+ 



x + a n + aj a n — cij 

2 \ f x + a, 
In 



In 



x + aj 



Or. 



v x + a n + aj aj — a n 

In this order of perturbation theory the same result is valid for the CP transformed processes. 
For the scattering I + I <-> h + H\ one has 



(E.14) 



E (^) 



x 



x + a,- x + 2a,- 



In 



x + a. 



+ 2^ Re 



j<n 



A^ 



71 J 



+ 



x + a n + a,- a n — a 5 



In 



x + a 



a, 



+ 



x + a n + a,- + a,- - a n 



In 



x + a r , 



a T , 



(E.15) 



The 2^3 process H 2 + q ] <-> Z t + Z t + Z7 C gives 



— - I 1- -In 



1 , / x + 2a, 



In 



x + a, 



a. 



+ 



£ Re [(At A,) 



^/ Oi n CLj 



x + 2a, 



5 Sp 



x + aj 



x + 2a 



>^±^+ln 

a n a, - 



x + a n + a, 



In 



x + a 



a. 



+ ( 2 ^±^- + In 



x + a n + a, 



In 



x + a r 



(In 



(E.16) 



+ Sp 



x + a n + a,- 



Sp 



x + aj 
x + a„ + a, 



+ Sp 



an 



x + a n + a. 



Sp 



x + an 
x + a n + a 3 j j 



For the related transition I + I ^ U c + q + H ] 2 we have 



ai 16) (x) 



3a« 
16vr 2 



^ \ /jj x [aj \ aj J 6 
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1 - 2 In 



x + 2a 7 - \ \ / x + a 1 
In 



(E.17) 



+ 2^Re 



aU,. 



nj 



^J(X n CLj 



' H " a 3 j^l ^ x + a n + ^ ^ ^ X + a j \ , G / X + CLj 



+ S P 



^±^ + ln( X + a " + a 0)ln(^±^)+S P (-^i) + ^ + iln 2 fe 
aj -a„ V a,- I I \ a n J \ a,- / 6 2 V a 



There are some 2^2 processes left which do not violate lepton number but simply transform 
leptons into scalar leptons, like in the process I + 7 <-> /i + H2 



-x , / x + a. 
+ In 



x + a,- 



+*E 



n,3 



In 



x + a 



+ °" In 



x + a r , 



or in the similar process I + He, *-* h + I 



(E.18) 



x + 2a,j . fx + a 7 - 
-2 H i In ' 



x 



+ 2 E 



n,3 





2 - 







-1 + 



Oj x + aj j ^ x + a,- \ , a n x + a n , ( x + a r 



+ 



In 



Finally, the last process gives 



(E.19) 



x + 2a,- / x + a, 
-2 H ^ In 



+*E 



aj x + aj ( x + a, \ a n x + a n ( x + a n 
-1 H — —In H In ' 



37 Cbj 



(E.20) 



E.2 Scattering off a Top or a Stop 

For the processes specified in fig. |<0 the reduced cross sections read 
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ar> = 3a u (Xt\ v 



(2) 



3a u (XlXu 



X — CLj 



jj x 



x — a» 



jj x 



2x — cij + 2cih x + 2a/j / x — aj + ah 
H In 



x — a 7 • + a/j x — aj 



a h 



x — aj , / x — aj + ah 
- 1 h In ' l 



x — a,; + 2a/j 



(3) 



3a u (aJ,A 



(7 



(i) 
u 



jj V £ 



jj x 



x — 2a,- + 2a/j a,- — 2a^ /x — a, • + 
H — in ' 



x — aj + ah x — aj 



j 



a h 



(E.22) 
(E.23) 
(E.24) 
(E.25) 



To regularize an infrared divergence in the t-channel diagrams we had to introduce a Higgs-mass 

2 



a h :- 



JL 

Mi 



(E.26) 



In the calculations we have used the value \i = 800 GeV. 



The analogous processes involving a scalar neutrino (cf. fig. 3.10) give similar contributions 

2 



a 



a 



(6) 



(5) _ 3a M 
2 

3a u 
3a u 
3a u 
3a u 



(7) 



*( 8 ) 



AtA, 



A+A, 



AtA, 



a; — gl 



jj V a; 



x — a,- 



x - aj + 2a h , { x - aj + a h 
-2 H In 



x — o, 



a h 



jj 



x — aj , ( x — aj + ah 
+ In 1 



x — a 7 - + 2a/j 



jj X X 

a 



jj x 



x — aj I x — aj + ah 

+ In ' 



x — a, + au 



aj, 



(E.27) 
(E.28) 
(E.29) 
(E.30) 
(E.31) 



E.3 Neutrino Pair Creation and Annihilation 



With the abbreviations 



v ij 



j ij 



X(x,ai,aj) 



x — \/a. 



X I \/ Qjj \ / CL'j 



In 



x — ai — aj + \/\ 



X (X^ \J 



the reduced cross sections for the right-handed neutrino pair creation read 



^=^{( A Hj( Ak 



2 Re 



j'. 



y/aiaj (ai + aj) _ 
x (x — ai — aj) 



(E.32) 
(E.33) 



n ( i 



(E.34) 
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<7 (2) 



1 4tt 



a 



(3) 



At A J (Xtx u , 



2 ^ Oi + Oj 

Mi 1 J-ii 



l «7 ■'-'ii 



- 2 Re 


(aJa) 2 . 







x — a; — a 



.7' 



2 r 



2 Re 



L 



(E.35) 
(E.36) 



«T (4) 



J ' 47T 



2 ^— Qj + aj 



X ' X 



- 2 Re 


(a>a„);_ 







a ; a 



(E.37) 



For the scalar neutrinos one has similarly 

1 



JV?iV c 

« J 



JV C JV C 



JV C JV9 



N C N C 



47T 



At A J.. AtA, 



2 r — | x-a i -a j T 



i{(AtA„)..(A»A„) j J^-2Re[(AtA„ 



a,- 



-Li 



X 



*7 [ ' 



1 

47T 

1 

4ir 





2 - 







" — V Aii + 

X 



3? G/^ T 

u n I Ljj 



Ay - 2 Re 



AtA,).. 



J'. 



yj 'CL{ CLj 



^3 



(E.38) 
(E.39) 
(E.40) 
(E.41) 



For the diagrams involving one neutrino and one sneutrino (cf. fig. 3.14) one finally has 



a 



(i) 



1 



N 3 N t Air 



- AU„ At A 



X + ai ~ a iL ij -2Re 

a x 



31 



^/OjOj x + a,{ — aj _ 

X X fl^ 



n ( ' 



(E.42) 



a (2 L = J_ 
^ 4vr 



aIaJ . 



J' 



x + di — aj 



Uj - 2 Re 



j i 



^Ja~iOj x + di — aj 

X X ttj 



Uj > ■ (E.43) 



Appendix F 



Reaction Densities 



In general the reaction densities corresponding to the reduced cross sections discussed in appendix |E] 
have to be calculated numerically. However, there exist some interesting limiting cases where one 
can calculate them analytically. Since thermal averaging of reduced cross sections via eq. ( D.16| ) in- 
volves modified Bessel functions, we start by summarizing a few useful formulae for Bessel functions 
before discussing the reaction densities. 



F.l Bessel functions 



Modified Bessel functions with different indices are related via recursion relations [53] 



xK v _i(x) + xK v+ i(x) = 2uK u (x) 

K v - 1 {x)-K v+1 (x)=2-fK v (x) 

ax 



For integer index Bessel functions have the following series representation 
2 



K„(x) 



1 n— 1 

EH 



fc=0 



k (n- k- 1)! 



+ 



Z \n+2k 



k=0 



k\(n + k)\ 



where ip denotes the derivative of the logarithm of the Gamma function 

ijj(x) = — — InT(x) . 
ax 



(F.l) 
(F.2) 



(F.3) 



(F.4) 
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For integer argument it reads 

n— 1 -. 

^( n ) = - 7B + £ (F.5) 
fc=i 

where 7 E = 0.577216 is Euler's constant. Hence, the leading terms of the series are given by 

K (x) = ln0) - 7 E + ... , (F.6) 
(n - IV /2\ n 

K„(z) = - +.-., forn>l. (F.7) 



2 V^, 

The asymptotic expansion of modified Bessel functions reads 

K (x) - ,/le- f -J— r (" + * + *) (F 8) 

fe=0 v ' 

i.e. to leading order all Bessel functions have the same asymptotic behaviour, 

K v {x) = J^e- + ... . (F.9) 

Furthermore, when evaluating reaction densities according to eq. ( p. 16 ), one has to compute 
integrals involving Bessel functions. In the following we compile the integrals that we have used in 
our calculations. 
Bessel functions and powers 

oo 

I dx x»K u (ax) = 2^a-^T (^f^) r {^fA (F.10) 
o 

for Re(l + n ± v) > and Re(a) > , 

l 

f dx x u+1 KJax) = 2 u a- v - 2 T(v + l) - -K u+1 (a) for Re(v) > -1 , (F.ll) 
J a 
o 

l 

f dxx 1 -"KJax) = 2- y a u - 2 T(l - u) - -K„_i(a) for Re(v) < 1 , (F.12) 
J a 
o 

oo 

dxK^zx) = -K Q (za) , (F.13) 

oo 

a 2 

dxx 2 K 1 (zx) = —K 2 (za) , (F.14) 
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oo 

f 2 2 

dxx Ki(zx) = — , 

o 

oo 

J dx xKq(zx) = -Ki(za) , 

a 
oo 

/dxxK (.zx) = -7T , 
z l 



oo 

y \/xKi(z\/i) = 2-K2(z\/a) , 

a 
oo 

y <ix \fx Ki (z\/x ) = — , 



oo 

/l 2 
y/x z 



(F.15) 



(F.16) 



(F.17) 



(F.18) 



(F.19) 



(F.20) 



Bessel functions and logarithms 



oo 

f 1 1 

/ dx x ln(x)Ko(zx) = — Ko(za) + a ln(a)-Ki(za) , 

J z z 



(F.21) 



oo 

y fix x ln(x)Ko(2x) = — 



-7e + In - 



(F.22) 



oo 

j dx^\n{^)K l {z^)= A + ^lnQ 



Kof^Va) 



+4^| [l + lng)]^^), 



oo 

y fife \/xln (j^j Ki(z\fx) 



(F.23) 
(F.24) 



In all these integrals we have always assumed that a, b and z are real and positive. 
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F.2 Lepton Number Violating Scatterings 



(i) 

In the Boltzmann equations we do not need every reaction density 7jy , i = 1, . . . , 19 separately 
(cf. sect. [O]). We only have to consider the combined reaction densities 





= 2 7 «+#- 




= 7l 3) +7l 4) - 




= 37if ) +7 r ) 




= 4 7 if ) + 27l 8) 



7i 4) +7i 6) +7i 7) +2 7 i 12) +7i 14) 



IN 



(7) , (14) 

7w + 7iV 



+7r+6 7 r 



,(io) 



„(9) 



y r + 2 7 r ) +7r ) +7r+67r ) 



(F.25) 
(F.26) 
(F.27) 
(F.28) 



For low temperatures, i.e. z S> , the dominant contribution to the integrand of the reaction 

densities comes from small centre of mass energies, i.e. x <C a,-. In this limit the reduced cross 
sections a$ for the (L + L) violating or conserving processes behave differently. For the (L + L) 
violating scatterings (i = 1, . . . , 5, 8, 10, 12, ... , 16) one finds 



(xx for x <C a,- , 



(F.29) 



while one has 

af oc x 2 for x <C a, (F.30) 

for the (L + L) conserving processes (i = 6, 7,9, 11, 17, 18, 19). In diagrams with an intermediate 
neutrino this different behaviour is due to the different chiral parts of the fermionic propagator con- 
tributing to the scatterings. (L + L) violating processes contain the chirality violating propagators 
y/aj/(x — a,j), whereas (L + L) conserving processes depend on the chirality conserving propagator 
y/x/(x — a,j). For diagrams with an intermediate scalar neutrino the corresponding kinematical 
factors originate in the couplings of sneutrinos to different initial and final states. 

Hence, the reaction densities can be calculated analytically in this limit and one finds 



Mf 1 

7T 5 



E 

3 



XtX u Y — + VRe 



19 



(F.31) 



Mf 1 

7T 5 Z 6 



A.yJ (l n (Lj 



(F.32) 



IC — _R _a 



E 



x 4 \ 2 1 / 27a u 



+ 



E 

j<n 



+ 



18a u 



7T 



(F.33) 
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ID 



7T 5 7T Z 6 



147 



nj'J 16^/a n aj 



(F.34) 



For high temperatures, i.e. z <C l/y/Oj , we can use the asymptotic expansions of the reduced cross 
sections to compute the reactions densities and we get 



Mf 1 
6471-5? 



47T 



n > 3 



a n — dj \ a 



+ ( 2 +^)E 



(F.35) 



IB 



Mf 1 

647T 5 "? 



— in 1 



«n ~ a j \ a 



j<n 



(F.36) 



7 Ai 



1 

32tt 5 F 



E 0^) 



45q m 9q m Ioj 
8vr + ~8~V~£. 



4 + 



27a M 
2vr 



In 



7e 



+E 

3<n 



2 + -2 a„W cijCj + a 

(a n -Oj) v v / 



+ 8 + 



36a u \ / a n 



In 



/ \ Q"n ~ a j 



+ 



In 



(F.37) 



- 7e 



i D 



Mf 1 
32vr 5 ? 



EM 



27a,, 39a,, fa7 ( , 27a. 

+ — ^ + — —,Hr+ 4 + 



+ 



E 



n,3 



nj 



4 - 8 7e + 



8vr 8 V Ci 



a„ In 



/ a n z 



2vr 



— cij In 



In 



z./o, 



7e 



(F.38) 



7T ^ 



sj a n aj 



7 1 



In 



+ 



57T 



2a n - aj \ aj J ( a „ - aj) 
i.e. these reaction densities are proportional to T 4 at high temperatures, as expected on purely 
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dimensional grounds. 

For intermediate temperatures z ~ 1/ \faj the reaction densities have to be computed numer- 
ically. This becomes increasingly difficult in the narrow width limit, where 1/Dj(x) has two very 
sharp peaks. However, in the limit Cj —* the two peaks in 1/Dj(x) cancel each other, since they 
have a different sign, while the peaks in 1/Dj(x) add up. Therefore, the terms proportional to 
1/Dj(x) or l/Dj(x)D n (x) with n ^ j can be neglected in the narrow width limit, while 1/Z)|(x) 
can be approximated by a 5-function 

1 7T 



D](x) ~ 2^0-75 

-2 



5(x-aj) . (F.39) 



An analogous relation holds for 1/Dj (x). 

These relations allow to calculate the contributions from the s-channel diagrams to the reaction 
densities analytically in the limit Cj — ► 0, while the contributions from the t-channel diagrams can 
easily be evaluated numerically. 



F.3 Stop Neutrino Scatterings 



The reaction densities 7 t for the interaction of a (s)neutrino with a top or a stop can also be 
calculated analytically in the limit of high temperatures z -C l/y777 . For the s-channel processes 
one finds 



7t 



(0) 



64vr 4 

„(o) 



K 2 (zy/aj) 



33 



It} = 2 %' > ~l 



(5) _ JO) 



(F.40) 
(F.41) 



For the t-channel reaction densities one has analogously 



7t 



(i) _ 3a u Mf 
8vr 4 



1 



33 f 



z 2 aj 



K [zy/aj) + 



ln[5i 

ah 



1 Kt (z^fal 



(2) 



3a u Mf 



8vr 4 



ln(^)K 2 (^) 



7t 



(4) 



27, 



(0) 



X 1 ) 



(7) 



(2) 



(F.42) 

(F.43) 
(F.44) 



7^ and j^' are several orders of magnitude smaller than the other 7^ for small z and can therefore 
be neglected at high temperatures. 



By using the series expansions (F\6) and ( |F.7| ), one sees that the processes with a higgsino in 



the i-channel, i.e. "ft , 7^ , 7^ and Jp' , behave like T 4 In (T/Mj) at high temperatures, whereas 
the other reaction densities are proportional to T 4 . 



S 7 ) 
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F.4 Pair Creation and Annihilation of Neutrinos 

In the Boltzmann equations we only need certain combinations of reaction densities which can 
easily be evaluated for high temperatures 



2 2 




i.e. these reaction densities are proportional to T 4 In (T/(M, + Mj)) at high temperatures. 



Bibliography 



[1] E. E. Cummings, 73 poems, (Faber and Faber, London, 1990) 

[2] For a discussion and references, see E. W. Kolb and M. S. Turner, The Early Universe, 
(Addison- Wesley, Redwood City, CA, 1990) 

[3] A. D. Sakharov, JETP Letters 5 (1967) 24 

[4] For a recent review, see R. H. Brandenberger, Particle Physics Aspects of Modern Cosmology, 
to be published in Field Theoretical Methods in Fundamental Physics, ed. by Choonkyu Lee, 
Mineumsa Co. Ltd., Seoul, 1997, Preprint BROWN-HET-1067, hep-ph/9701276 

[5] L. Kofman, A. Linde and A. A. Starobinsky, Phys. Rev. D 56 (1997) 3258 
I. Zlatev, G. Huey and P. Steinhardt, Phys. Rev. D 57 (1998) 2152 

[6] I. Affleck and M. Dine, Nucl. Phys. B 249 (1985) 361 

[7] see, e.g., L. Randall, in Electroweak Interactions and Unified Theories, Proc. of the XXXth 
Rencontres de Moriond, ed. J. Tran Thanh Van, (Editions Frontieres, 1995) 

[8] For a review and references, see V. A. Rubakov and M. E. Shaposhnikov, Phys. Usp. 39 (1996) 
461; 

[9] For a discussion and references, see K. Jansen, Nucl. Phys. B (Proc. Supp.) 47 (1996) 196 

[10] V. A. Kuzmin, V. A. Rubakov and M. E. Shaposhnikov, Phys. Lett. B 155 (1985) 36 

[11] T. Yanagida, in: Workshop on unified Theories, KEK report 79-18 (1979) p. 95; 

M. Gell-Mann et al., in: Supergravity, Proc. of the Stony Brook Workshop 1979, eds. P. van 
Nieuwenhuizen and D. Freedman (North Holland, Amsterdam, 1979), p. 315 

[12] M. Fukugita and T. Yanagida, Phys. Lett. B 174 (1986) 45 

99 



100 



BIBLIOGRAPHY 



[13] M. A. Luty, Phys. Rev. D 45 (1992) 455 

[14] M. Plumacher, Z. Phys. C 74 (1997) 549 

[15] W. Buchmiiller and M. Plumacher, Phys. Lett. B 389 (1996) 73 

[16] B. A. Campbell, S. Davidson and K. A. Olive, Nucl. Phys. B 399 (1993) 111 

[17] M. Plumacher, hep-ph/9704231, to appear in Nucl. Phys. B 

[18] L. Covi, E. Roulet and F. Vissani, Phys. Lett. B 384 (1996) 169 

[19] A. Yu. Ignatev,V. A. Kuzmin and M. E. Shaposhnikov, JETP Lett. 30 (1979) 688; 
F. J. Botella and J. Roldan, Phys. Rev. D 44 (1991) 966; 

J. Liu and G. Segre, Phys. Rev. D 48 (1993) 4609; Phys. Rev. D 49 (1994) 1342 

[20] M. Flanz, E. A. Paschos and U. Sarkar, Phys. Lett. B 345 (1995) 248; ibid. B 384 (1996) 487 
(E); 

L. Covi, N. Rius, E. Roulet and F. Vissani, Phys. Rev. D 57 (1998) 93 

[21] M. Flanz, E. A. Paschos, U. Sarkar and J. Weiss, Phys. Lett. B 389 (1996) 693; 
L. Covi and E. Roulet, Phys. Lett. B 399 (1997) 113; 
A. Pilaftsis, Nucl. Phys. B 504 (1997) 61; Phys. Rev. D 56 (1997) 5431 

[22] W. Buchmiiller and M. Plumacher, Phys. Lett. B 431 (1998) 354 

[23] M. Veltman, Physica 29 (1963) 186 

[24] E. Roulet, L. Covi and F. Vissani, hep-ph/9712468 

[25] S. J. Gates, M. T. Grisaru, M. Rocek and W. Siegel, Superspace or One Thousand and One 
Lessons in Super symmetry, (Benjamin/Cummings, Reading, MA, 1983) 

[26] H. P. Miles, Phys. Rep. 110 (1984) 1 

[27] M. Sohnius, Phys. Rep. 128 (1985) 39 

[28] A. Salam and J. Strathdee, Nucl. Phys. B 76 (1974) 477 

[29] P. P. Srivastava, Super symmetry, Superfields and Supergravity, (Adam Hilger, Bristol, 1986) 

[30] J. Wess and J. Bagger, Supersymmetry and Supergravity, (Princeton University Press, Prince- 
ton, NJ, 1992) 



BIBLIOGRAPHY 



101 



[31] B. Zumino, Nucl. Phys. B 89 (1975) 535 

[32] M. T. Grisaru, W. Siegel and M. Rocek, Nucl. Phys. B 159 (1979) 429 

[33] J. Iliopoulos and B. Zumino, Nucl. Phys. B 76 (1974) 310 

[34] E. Witten, Nucl. Phys. B 188 (1981) 513 

[35] S. Dimopoulos and L. Susskind, Phys. Rev. D 18 (1978) 4500 

[36] A. D. Dolgov and Ya. B. Zeldovich, Rev. Mod. Phys. 53 (1981) 1; 

E. W. Kolb and S. Wolfram, Nucl. Phys. B 172 (1980) 224, ibid. B 195 (1982) 542 (E) 

[37] W.-Y. Keung and L. Littenberg, Phys. Rev. D 28 (1983) 1067 

[38] S. Yu. Khlebnikov and M. E. Shaposhnikov Nucl. Phys. B 308 (1988) 885 

[39] L. E. Ibahez and F. Quevedo, Phys. Lett. B 283 (1992) 261 

[40] W. Buchmiiller and D. Wyler, Phys. Lett. B 249 (1990) 458 

[41] L. Wolfenstein, Phys. Rev. Lett. 51 (1983) 1945 

[42] W. C. Haxton, Annu. Rev. Astron. Astrophys., 33 (1995) 459 

[43] L. Wolfenstein, Phys. Rev. D 17 (1978) 2369; 

S. P. Mikheyev and A. Y. Smirnov, Nuovo Cim. 9C (1986) 17 

[44] E. Witten, Phys. Lett. B 91 (1980) 81; 

Z. Berezhiani, in: High Energy Physics and Cosmology, Proc. of the 1995 ICTP Summer 
School, Trieste, eds. E. Gava, A. Masiero and K. S. Narain (World Scientific, Singapore, 
1997), p. 618 

[45] T. Moroi, Ph. D. Thesis, hep-ph/9503210 

[46] Super-Kamiokande Collaboration, Y. Fukuda et al., hep-ex/9803006; 

E. Kearns, Talk at News about SNUS, ITP Workshop, Santa Barbara, Dec. 1997, 
http : //www. itp .ucsb . edu/online/ snu/kearns/ 

[47] H. Murayama and T. Yanagida, Phys. Lett. B 322 (1994) 349; 

M. Dine, L. Randall and S. Thomas, Nucl. Phys. B 458 (1996) 291; 

G. Lazarides, R. K. Schaefer and Q. Shafi, Phys. Rev. D 56 (1997) 1324 



102 



BIBLIOGRAPHY 



[48] B. A. Campbell, S. Davidson, J. Ellis and K. A. Olive, Phys. Lett. B 297 (1992) 118; 
J. M. Cline, K. Kainulainen and K. A. Olive, Phys. Rev. D 49 (1994) 6394 

[49] D. Boyanovsky, I. D. Lawrie and D.-S. Lee, Phys. Rev. D 54 (1996) 4013; 
A. Riotto, hep-ph/9712221; 

I. Joichi, Sh. Matsumoto and M. Yoshimura, hep-ph/9803201 
[50] G. Passarino and M. Veltman, Nucl. Phys. B 160 (1979) 151 
[51] A. Denner, Fortschr. Phys. 41 (1993) 307 

[52] A. Denner, H. Eck, O. Hahn and J. Kiiblbeck, Phys. Lett. B 291 (1992) 278; Nucl. Phys. B 
387 (1992) 467 

[53] I. S. Gradstein and I. M. Ryshik, Tables of Series, Products and Integrals, (Harri Deutsch, 
Thun, 1981); 

A. Erdelyi, W. Magnus, F. Oberhettinger and F. G. Tricomi, Higher Transcendental Functions, 
vol. II, (McGraw-Hill Book Company, 1953) 



103 



Acknowledgments 

I would like to thank my advisor, W. Buchmiiller, who suggested this investigation, for continuous 
support and a close and fruitful collaboration. 

I would also like to acknowledge highly instructive discussions with W. Beenakker, J. Hein, 
G. Hiller, A. Jakovac, S. Moch, T. Plehn and U. Nierste. 

Finally, I would like to thank all my colleagues and friends at DESY for the enjoyable working 
atmosphere. 



